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We perform a QCD analysis of the exclusive production of two mesons in 7*7 collisions in 
the kinematical domain of large photon virtuality Q and small hadronic invariant mass W . This 
reaction is dominated by a scale invariant mechanism which factorizes into a pcrturbative subprocess, 
7*7 — * QQ or 7*7 — * 99, an< A a generalized two-meson distribution amplitude. We develop in detail the 
phenomenology of this process at e + e~ colliders. Using a simple model for the two-pion distribution 
C**) , amplitude, based on its general properties, we estimate the cross section for the kinematics accessible 

at BABAR, BELLE, CLEO and LEP. 
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Exclusive hadron production in two-photon collisions provides a tool to study a variety of fundamental aspects 
of QCD and has long been a subject of great interest (cf., e.g., J|-||] and references therein). Recently a new facet 
r-j • of this has been pointed out, namely the physics of the process 7*7 —> 7T7t in the region where Q 2 is large but W 2 
O |' small j|. This process factorizes |||| into a perturbatively calculable, short-distance dominated scattering 7*7 — > qq 
or 7*7 — > gg, and non-perturbative matrix elements measuring the transitions qq — > 7T7r and gg — > tttt. We have 
called these matrix elements generalized distribution amplitudes (GDAs) to emphasize their close connection to the 
distribution amplitudes introduced many years ago in the QCD description of exclusive hard processes 0. 

Indeed it is instructive to consider 7*7 — ► 7T7t as a generalization of the process 7*7 — ► 7r°, where the distribution 
amplitude of a single pion appears. The 7-7r transition form factor has been the subject of detailed theoretical 
studies |@j. The experimental data |J are well reproduced by a description based on QCD factorization and provide 
S— i ' one of the best constraints so far on the form of the single-pion distribution amplitude. 

From a different point of view 7*7 — > 7T7t is the crossed channel of virtual Compton scattering on a pion. The 
kinematical region we consider here is closely related to deeply virtual Compton scattering (DVCS), which has 
attracted considerable attention in the context of skewed parton distributions pQ] . 

Our reaction can also be seen as the exclusive limit of a hadronization process. The hadronization of a 55-pair 
originating from a hard, short-distance process such as a 7*7 collision is usually formulated in terms of fragmenta- 
tion functions which describe in a universal way semi-inclusive reactions, specifically the transition from a quark or 
antiquark to a final-state hadron when one integrates over all final states containing this hadron. We specialize here 
to the case where the final state consists of two mesons with specified four-momenta, and nothing else. 
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I. INTRODUCTION 
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Like other hadronic matrix elements the GDAs are process independent. It has recently been pointed out that 
they occur in the hard exclusive process j*p — > irirp, where the pion pair is or is not the decay product of a p meson, 
and that the analysis of that reaction would benefit from the measurement of the two-pion GDA in 7*7 — > tttt. 

All these aspects lead us to consider GDAs as a promising new tool for hadronic physics, which may be used to 
unveil some of the mysteries of hadronization and the confining regime of QCD. The process 7*7 —> tttt is well suited to 
access these quantities experimentally. In the present paper, we develop in detail the phenomenology of this reaction 
and emphasize the feasibility of its investigation at existing e + e~ colliders. 

In Sect. we discuss the kinematics of our process, recall its main properties in the factorization regime we 
are interested in, and sketch the crossing relation between 7*7 — ► tttt and deep virtual Compton scattering. In 
Sect. |n| we list the general properties of generalized distribution amplitudes and in particular review their QCD 
evolution equations. These properties lead us to construct a simple model of the two-pion GDA, which is described 
111 Sect. [Ev[ Section [v| gives a comparison between one-pion and two-pion production in 7*7 collisions. Relations 



with the inclusive production of hadrons, commonly described by the photon structure function, are discussed in 



Sect. VI. The phenomenology of our process in e7 collisions is described in detail in Sect. VII, with special emphasis 



on the information contained in angular distributions and in the interference with the bremsstrahlung mechanism. In 



Sect. VIII we give estimates for the cross section for various experimental setups at existing e + e~ colliders. Section IX 
contains our conclusions. In Appendix [a| we specify our sign conventions for pion states, and in Appendix |b| we discuss 
what additional information can be obtained with polarized beams. 

II. THE PROCESS 7*7 -> tttt 
A. Kinematics in the 7*7 center of mass. 

The reaction we are interested in is 

e(fc)+7(</) -» e (fc')+7r , '(p)+7r'V). (!) 
where four-momenta are indicated in parentheses. We further use 

q = k-k', Q 2 = -q 2 , P = p + p', W 2 = P 2 . (2) 

The pions may be charged (i = +, j = — ) or neutral (i = j = 0), and the lepton e may be an electron or a positron. 
Scattered with large momentum transfer this lepton radiates a virtual photon j*(q), and for the 7*7 subprocess we 
introduce the Bjorken variable 

- Q " (3) 



2q -q> Q 2 + W 2 ' 

In e + e _ collisions the photon 7(9') can be obtained by bremsstrahlung from the other beam lepton, so that the 
overall process is 

e(k) + e(i) -» e(fc') + e(Z') + tt>) + n j (p') (4) 

with q' = 1 — V . In the spirit of the equivalent photon approximation we approximate q' 2 as zero and the momenta q' 
and I as collinear. We write E\ — k°, E2 = 1° and q'° = X2 1° for the energies in the laboratory frame.fj] For the cm. 
energies of the ee and ej collisions we have 

s ee = (k + l) 2 , s el = (k + q') 2 = x 2 s ee . (5) 

Let us now discuss the kinematics in the 7*7 center of mass frame. We use a coordinate system with the z axis 
along q, and with x and y axes such that p lies in the x-z plane and has a positive x component, i.e., 

g=(g ,0,0,|q|), p= (p°,|p|sin#,0,|p|cos#), (6) 



1 We neglect the small finite crossing angle between the beams at BELLE, so that in our parlance the lepton beams are 
collinear in the "laboratory frame". 
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FIG. 1. The kinematics of e(k) + 7(9') — ► e(fc') + vr'(p) + 7r J (p') in the center of mass of the pion pair. 

where we have introduced the polar angle 9 of p. Another natural variable for our process in this frame is the azimuth 
ip of k', which is the angle between the leptonic and the hadronic planes, cf. Fig. [y. In terms of Lorentz invariants 
these angles can be obtained from 

cosfl — 



COS if 



(3{Q 2 + W 2 Y 

2k ■ (p' - p) (Q 2 + W 2 ) + (3cos6 [Q 2 {s ei -Q 2 -W 2 )-s ei W 2 ] 



2/3sin6i v /se 7 Q 2 W 2 (s ei -Q 2 -W 



4e af 3 jS (p + p') a p Wq 5 

sm<p= . = (7) 

P SnOy/Sery Q 2 W 2 { Se ~ t - Q 2 - W 2 ) 



with eoi23 = +1 and the velocity 



Am 2 
~W~ 2 



/3 = \/l-Tf7f W 



of the pions. A further quantity we will use is the usual y-variable for the e'y collision, 



y = W = Q 2 + W 2 



which can be traded for 



l-y + 2/72' 



(10) 



the ratio of longitudinal to transverse polarization of the virtual photon 7* (q) . 

We finally define light cone components — (a ± a 3 )/V2 for any four- vector a and introduce the fraction 

c _ P + _ l + /3cos0 
of light cone momentum carried by 7r l (j>) with respect to the pion pair. 

B. Factorization at large Q 2 and small W 2 

Let us briefly review how 7*7 — > rnr factorizes in the kinematical regime we are interested in. Firstly, we require 
Q 2 to be large compared with the scale A 2 ~ 1 GeV 2 of soft interactions, thus providing a hard scale for the process. 
Secondly, we ask W 2 to be small compared with this large scale Q 2 . In this regime the dynamics of the process is 
conveniently represented in the Breit frame, obtained by boosting from the 7*7 center of mass along the z axis. The 
spacetime cartoon of the process one can derive from power counting and factorization arguments is shown in Fig. ||. 
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In the Breit frame the real photon moves fast in the negative z direction and is scattered into an energetic hadronic 
system moving in the positive z direction. The hard part of this process takes place at the level of elementary 
constituents, and the minimal number of quarks and gluons compatible with conservation laws (color etc.) are 
produced. At Born level one simply has 7*7 — > qq, but through a quark box the photons can also couple to two gluons. 
Each quark or gluon carries a fraction z or 1 — z of the large light-cone momentum component P + . Subsequently the 
soft part of the reaction, i.e., hadronization into a pion pair, takes place. 

At leading order in as the amplitude is given by the diagram of Fig. I (a) and the one where the two photon 
vertices are interchanged. One calculates for the hadronic tensor Q 

TV = ij d±x e-^ (7rCp)7rCp')| TJgJx) J^(0) |0) = -tf £ | j\z i£L $-( 2j C , W 2 ), (12) 

where gt^ denotes the metric tensor in transverse space (g 11 = — 1). The sum on the r.h.s. runs over all quarks flavors, 
e q is the charge of quark q in units of the positron charge e, and eJ^ m (x) is the electromagnetic current. While the 
expression of the hard subprocess 7*7 — > qq is explicit in Eq. (|l2|), the soft part of 7*7 — > tttt is parameterized by the 
generalized distribution amplitude 

*™{z, C, W 2 ) = J ^- e -^r + x-) (n{p)n{p')\ q(x-)>y + q(0) |0) (13) 

for each quark flavor q. We work in light cone gauge A + = 0, otherwise the usual path ordered exponential of gluon 
potentials appears between the quark fields. $ 9 depends on the light-cone fraction z of the quark with respect to the 
pion pair, on the kinematical variables £ and W 2 of the pions, and on a factorization scale. The latter dependence, 
not displayed in Eq. (|l3|), will be discussed in Sect. IIIB. 



In Eq. (|12|) a scaling behavior for our process is manifest: at fixed £ and W 2 the 7*7 amplitude is independent of 
Q 2 , up to logarithmic scaling violations from radiative corrections to the hard scattering and from the evolution of 
the two-pion distribution amplitude. This scaling property is central to all processes where a factorization theorem 
holds, and it is the basic signature one looks for when testing whether the asymptotic analysis developed here applies 
to an experimental situation at finite Q 2 . There will of course be power corrections in A/Q and W/Q to this leading 
mechanism. Examples are the hadronic component of the real photon, and the effect in the hard scattering of the 
transverse momentum of the produced parton pair. We note that the crossed channel, i.e., virtual Compton scattering 
has been analyzed in detail within the operator product expansion |p|, p^|Jl3| , which may provide a framework for a 
systematic study of higher twist effects. 

Contracting the hadronic tensor ( fl2"| ) with the photon polarization vectors we see that in order to give a nonzero 
7*7 — > tttt amplitude the virtual photon must have the same helicity as the real one. As in the case of deep virtual 



Compton scattering this is a direct consequence of chiral invariance in the collinear hard-scattering process [|14 15 
and is valid at all orders in as- In the case of the 7*7 — > gg subprocess the photon helicities can also be opposite 
p6[ . In any case the virtual photon must be transverse. As a consequence nonleading twist effects can be studied in 
the amplitude for longitudinal 7* polarization, without any "background" from leading twist pieces. We will develop 



in Sect. VII how the different 7*7 helicity amplitudes are experimentally accessible. 
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(a) (b) 




FIG. 3. (a) Factorization of the process 7*7 — > tvtv in the region Q 2 3> W 2 ,A 2 . The hard scattering is shown at Born 
level, with a second diagram being obtained by interchanging the two photon vertices. The blob denotes the two-pion GDA. 
(b) Crossing relates this process to deep virtual Compton scattering, 7*7r — > 771-. The blob now denotes the (skewed) quark 
distribution in the pion. 



As we already mentioned, there is a close analogy of two-pion production in the region Q 2 ^> W 2 , A 2 with the one- 
pion channel, commonly described in terms of the 7-71 transition form factor. There again a factorization theorem 
holds, which allows the hadronic tensor T^ u to be expressed in terms of the single-pion distribution amplitude (j)^ as 

T^ = ij d*x e-** (tt°| TJ? m (x)JZ m (0) |0) = £ | j\z (14) 

to leading order in as, where et^ is the antisymmetric tensor in transverse space (elf = 1) and 

r q {z)=i J ^-e~" ( p+ ^)(^(P)|g>-)7 + 75 9(0)|0). (15) 

Notice the different Dirac structures in the matrix elements ([l3|) and (|l5|), due to the different parity transformation 
properties of one- and two-pion states El . 

The theoretical analysis of this process has been highly developed [0 . Its generalization to the production of rj and 
t]' is also important, in particular with respect to the SU(3) flavor structure of the QCD evolution equations and the 
mixing of the quark singlet and gluon channels [p"7| . In Sect. ^ we will further compare the production of a single 
pion with that of a pion pair. 



C. Relation with deep virtual Compton scattering and parton distributions in the pion 



The process 7*7 — > tttt at large Q 2 and s <C Q 2 is related by s—t crossing to deep virtual Compton scattering on 
a pion, i.e., to 7*7r — * jtt at large Q 2 and —t <C Q 2 ■ It turns out that factorization works in completely analogous 
ways for both cases, as is shown in Fig. |^. The non-perturbative matrix elements occurring in the Compton process 
are skewed parton distributions ]1C| ], defined in the pion case as |18 



c(P+z~) 



(16) 



with P = (p + p')/2. They have been recognized as objects of considerable interest and have triggered intensive 
theoretical and experimental work. The processes 7*7 — > ttt: and 7*7r — > 771 share many common features, from their 
scaling behavior and th e det ails of their helicity selection rules to the possibilities of phenomenological analysis, which 
we will develop in Sect. VII . 

The imaginary part of the forward virtual Compton amplitude, j*ir — > 7*7r, obtained from Fig. || (b) by replacing 
the 7 with a second 7*, gives the cross section for inclusive deep inelastic scattering, 7*7r — > X, where the ordinary 
parton distributions in a pion occur. 

As observed in Jl^ j it is useful to implement crossing at the level of moments in momentum fractions (z and £ for 
GDAs, x and £ for SPDs), which depend only on a factorization scale and a Lorentz invariant (s for GDAs, t for 
SPDs). The moments of GDAs and of SPDs are connected by analytic continuation in that invariant. In particular, 
analytic continuation to the point t — leads to mo men ts of the ordinary parton distributions in the pion, which we 
will use as an input for our model of GDAs in Sect. IV. 
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III. GENERAL PROPERTIES OF GDAS 



A. Charge conjugation and isospin properties 



Let us start by compiling some symmetry properties which will be useful in the following. For the quark GDAs 
( |l3| ) the invariance of strong interactions under charge conjugation C implies 

C, W 2 ) = -$7(1 - z, 1 - C, W 2 ). (17) 

It is useful to project GDAs for charged pions on eigenstates of C parity, 

$±(z, C, W 2 ) = \ (*, C, W 2 ) ± (z, 1 - C, W 2 )) , (18) 

so that 

*J + w ~ = *t( z > C> + C> v^ 2 ). (19) 

In the C even sector Eq. ( |l7j ) reduces to 

$+(z,C,VF 2 ) = -$+(l-z,C,^ 2 ). (20) 

Our process is only sensitive to the C even part of Qg +7r since the initial state two-photon state has positive C parity. 

Of course a 7r°7r° pair has positive C parity as well, so that has no C-odd part at all. 

Let us now turn to isospin symmetry. The C odd component of a two-pion state has total isospin 1=1, whereas 
its C even component contains both 1 = and 1 = 2 pieces. The quark operator in <i>J 7r has only components with 
isospin / = or I = 1. Hence it is a consequence of the leading twist production mechanism and of isospin invariance 
that in our process the pion pair is in a state of zero isospin, i.e., that no component with I = 2 is produced. Another 
consequence of isospin invariance is that 

t>f*°(z,(,W 2 ) = <l>+(z,(,W 2 ), (21) 

so that the production amplitudes for neutral and charged pion pairs are equal. Deviations from isospin symmetry in 
the present reaction would be interesting, but since one can expect them to be small we will assume isospin invariance 
to hold throughout the rest of our study. Isospin invariance also implies that 

= *5, *z = ( 22 ) 

so that in the C even sector we only need to know the SU{2) flavor singlet combination + <i>J. 
The connection between the notation 'I*^ -0 ' 1 of Polyakov |l9| and ours is 

*f° = 3+, <J>H =1 = <iV (23) 

We remark that the second term in Eq. (2.6) of Ref. jl9| should come with a minus sign |2(J. Our relation $j^ =1 = $« 
takes this correction into account. 

Notice that the signs in Eqs. ([H]) and (|2^) depend on the choice of relative phases in the definition of charged pion 
states. We specify our convention in Appendix |X| 



B. Evolution equation 



In the process of factorization generalized distribution amplitudes acquire a scale dependence in the same way as 
usual distributions do. This scale dependence can be computed within perturbative QCD, and there is nothing special 
with multiparticle states since the scale dependence is a property of the nonlocal product of fields under consideration, 
rather than one of a particular hadronic matrix element (see |21| for an approach exploiting this feature). The scale 
dependence of GDAs can be cast in the form of an ERBL evolution equation p2[ ], and the only complication in the 
channel we are concerned with here is the mixing of quark and gluon distribution amplitudes. The leading-logarithmic 
form of the evolution equations has been studied in detail for the parity-odd sector p7[ , where the relevant quark 
operator is ?7 + 7s <?• Our application to pion pairs leads us to consider the parity-even sector, where the quark 



G 



operator is q~f + q instead, see our remark after Eq. (|l5|). For completeness we give here the basic steps for deriving 
and solving the evolution equation in this channel, following the procedure outlined in |p3| . Taking into account the 
different normalization conventions we find agreement with the results of Baier and Grozin |p3 , who reported a sign 
discrepancy with Chase p5| for the gluon evolution kernel. 

We are then concerned with the generalized quark and gluon distribution amplitudes in A + = gauge: 

<S> q (z,(,W 2 ) = l^-e-" (p+ ^(n(p)n(p')\q(x-h + q(0)\0), 

$ 3 (z,C,VK 2 ) = ±- l^-e-^ p ^^(n(p)n(p')\F+^x-)F, + (0)\0), 

= z(l^z)P+ l^-e-^ p+x ^(^PMp'U^x-)A^O)\0). (24) 

Our gluon distribution amplitude $ g (z, £, W 2 ) coincides with $ G (z, £, W 2 ) introduced in From the definition (gj) 
one readily obtains 

$ g (z,(,W 2 ) = $ g (l-z,{,W 2 ), (25) 

and from C invariance one has 

$ s (z,C,T^ 2 ) = <i> 3 (l-z,l-C,^ 2 )- (26) 

Here we have given definitions for a two-pion state, but as stated above the evolution equation for DAs and GDAs is 
not specific to the details of the hadronic system. The considerations of this and the following subsection thus apply 
to any state in the parity even sector which has four-momentum P and total angular momentum J z — along the 
axis defining the light cone variables. 

We now study the evolution of the distributions for gluons and of quarks in the singlet combination of rif flavors. 
For convenience we introduce 

n f 

z 2 z 2 f G (z) = <f> g (z), (28) 

where we use the notation z = 1 — z. In the end we will return to the amplitudes $ g and $ s . 
The scale dependence is controlled by the parameter 

,2 ,.2n 2 > f <*s(l4) 



^■sl^J' (29) 

where as is the one-loop running coupling and (3\ = 11 — 2n//3. This parameter describes how the distribution 
amplitude evolves when one changes the factorization point from po to p. The evolution equation takes the form 



where / is a two-component vector 



9 f(z,0 = V*f = [ duV(z,u)f(u,0- (30) 



/-(£]. 



and V is the 2x2 matrix kernel 

— . 

\ Vgq Vgg 

To obtain the leading logarithmic evolution equation it is sufficient to consider one-loop corrections to the scattering 
amplitude. The latter is depicted in Fig. [| and has the form H * /, where H — (Hq, Hq) denotes the hard-scattering 
kernels. It turns out that in light cone gauge A + = the relevant one-loop diagrams consist of an insertion between 
H and / of the graphs shown in Fig. (a) to (e), supplemented by (renormalized) self-energy insertions on each 
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z z 

FIG. 4. The scattering amplitude H * f with / denoting the soft matrix elements and H the hard scattering kernels. 



It 




Z U 

a) 




c) 





b) 



d) 



e) 



FIG. 5. One-loop insertions, to be supplemented by self-energy insertions on every line appearing in Fig. |4|. The sum of all 
insertions gives the evolution kernel £V. We remark that the one-loop graph (e) must be multiplied by 1/2 to avoid double 
counting, u and z denote light cone plus momentum fractions, and k the loop four-momentum. 



line connecting H to / in Fig. [|. Calling the sum of these insertions the one-loop diagrams have the structure 
H*iV*f. 

The evolution from zeroth to first order of the generalized distribution amplitude may thus be written as 

/« (z) = /<°> (z) + £ ( du V(z, u) / (0) («). (33) 
Jo 

In the computation of the diagrams, the k~ integral is performed by the Cauchy method of contour integration in 
the complex plane, and £ is the result of the integral over transverse momentum from kt = Ho to kt = 



(34) 



Despite the presence of as in Eq. (p4|), ^ is not small if fi 2 ^> fig, and this signals the necessity of an all-order analysis. 
This analysis leads to the evolution equation, with the feature that V is the same matrix in Eqs. Q30| ) and (|33|). We 
refer the reader to the literature for a general discussion j|(J . 

The integration over k + may be reexpressed as an integral over the incoming light cone fraction u. The evolution 
kernels contain the remaining part of the dynamics, in particular they describe the change of light cone fractions from 
u to z. We get 



z, u) = Cf 



6(z-u)- 1 + 



1 



+ { 



u <-> u, z <-> z 



Vqg(z, it) = 2n/Tp 6(z — u) — (2z — u) — {u <-> u, z <-> z} 



Vgq(z,u) 
V G g(z,u) 



Cf 

zz 
C A 



(z — u) — {z — 2u) — {u «-> u, z <-> z} 

2 



2 2 



[(2z-l) 2 + (2u-l) 2 ] +{u 



<-> u, z <-> 



■n/Tp5(u — z), (35) 
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where the color factors are Cf — 4/3, Tp — 1/2 and Ca = 3. The subscript + stands for the + distributions, whose 
action on a function / may be expressed symbolically as 



V ' 1 + 



/(«) = [•••] (/(«)-/(*))• 



(36) 



The kernels ( p5| ) give the finite parts that remain after the cancellation of infrared divergences between graph (a), 
resp. (d), and quark self-energy, resp. gluon self-energy insertions. A simple way to obtain self-energy corrections is 
to notice their relation to parton splitting |23j , that is 



/£>(*) 





1 


-«/ 


SE 








1 


-</ 


SE 





i- if dxP GQ ( x ) 



l-£ J dx (^Pgg(x) + n f P QG (x] 
with the unregularized DGLAP splitting functions 



/£>(*), 



(37) 



P QG (:r)=T F [ a ; 2 + (l-.x) 2 ] 
P GG {x) = 2C A 



x 
x 



l-x 



1 — X 



+ x(l - x) 



(38) 



The integrals ( p7| ) are not defined in the limit a; — * 0,1, which is a manifestation of the infrared divergence of self-energy 
graphs. 



C. Solution 



We will now solve the evolution equation (f30|). Given our application we restrict ourselves to the C even parts 
of the quark distributions, the gluon distribution being of course even under C from the start. 
We look for solutions of the form 



f(z,0 = f(z)e-^. 



(39) 



To this end it is convenient to change variables, introducing y = 2u — 1 and x = 2z — 1, and to study the convolution 
of the matrix kernel V with 









„n-l I i 



where n is an odd integer to accommodate the symmetry properties (M) and feq). One finds 



Vqq * y n = -l QQ (n) x n + 0(x n ~ 2 ), V QG * y n ~ x 

Vgq * y" = ~lG Q {n) x n - x + 0(x n ' 3 ), Vqq * y"- 1 
with anomalous dimensions 

n+l 

k 



-7QG(n)x n + 0(x n - 2 ), 
-^G G (n)x n - 1 + 0(x n - 3 ), 



lQQ{n)=CF [\- {n + l)(n + 2) 



n+l \ 



fe=2 



7Qe(n) 
iGQin) 



-n f T F 
-2C F 



n 2 + 'in + 4 
n(n + l)(ra + 2)' 
n 2 + 3n + 4 



(n + l)(n + 2)(n + 3)' 



1 



lGG{n) = CU - 



6 n(n + l) (n + 2)(n + 3) 



fe=2 / 



(40) 



(41) 



(42) 
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Since for a given no the space of solutions with n < no is stable under the application of the kernel one can find 
polynomials p n (x) and q n —x(x) satisfying 



(43) 



* Pn = -lQQ{n)p n , V QG * = -jQc(n)p ni 

Vgq *p n = -jGQ(n) q n -i, Vgg *q n -i = -7GG(n)q n -i- 

The symmetry properties of the kernels 

(l-x 2 )V QQ (x,y) = (l-y 2 )V QQ (y,x), 
2C F (1 - x 2 ) Vqg(x, y) = n f T F {\ - y 2 ) 2 Vgq (y, x), 
(1 - x 2 ) 2 V GG (x, y) = (l- y 2 ) 2 V GG (y, x), 



(44) 



then imply that the (p n ) are orthogonal polynomials on the interval [—1, 1] with weight 1— x 2 , i.e., they are proportional 
to the Gegenbauer polynomials Cn^ 2 \x), whereas the (q n -i) are orthogonal on [—1, 1] with weight (1 — a; 2 ) 2 , that 
is, proportional to the Gegenbauer polynomials C n _i (x). To complete the identification it is necessary to take into 

r 3 / 2 ^ (5/21 

account the standard normalization of Gegenbauer polynomials. One finds that p n = C n and q n -i = C n L\ fulfill 
Eq. (|43]), provided one makes the replacements 



n 3 
7qgW -» 1qg(p) = g 7qg(")' 7gq(") -» 7gqM = ~7gqW- 



(45) 



The final step is to diagonalize the 2x2 anomalous dimension matrices for each value of n. The eigenvalues are 



r (±) 



1 



7qq» + 7gg0) ± \J[lQQ( n ) - 7gg(")] 2 + 4 7qg( 71 )7gq(«) 



and the eigenvectors of the kernel matrix read 

vit\x) 

where 



Ck 3/2 \x) 



fJ 



(±) 



r (±) 



(») 



7qg( 7 



The general C even solution of Eq. ( |30| ) may then be written as 

/M) = £ {A^\x)e-^ + AiM-Hx)e-^} 



(46) 



(47) 



(48) 



(49) 



odd n 



with integration constants A n ■ 

We now return to the amplitudes <fr q , $ g and explicitly express £ in terms of /i and «o- The key result of this 
section then reads 



J>+(z, U 2 ) = z(l-z) J2 A n (f, 2 )C^ 2 \2z-l), 

odd n 

$ s (z, M 2 ) = ^(l-^) 2 £ ^(M 2 )C^ /2) (2z-l) 



(50) 



with 



yi„(" 2 )-4 +) 



«s(» 2 ) 



odd n 



«s(M 2 ) 



if'-' 



(51) 
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and exponents = 2T^/f3i, which are positive except for k[ ^ — 0. For n/ = 2,3,4, one explicitly finds 



A"} 



( +) _ 32 + 6nj 



0.51, 0.62, 0.75, 



K. 



(-) 



0.71, 0.76, 0.82, 



K 



(+) 



1.45, 1.64, 1.85. 



(52) 



99- 6n f 

^From Eq. (|5l]) we easily see that the integration constants A^' depend on the starting scale /j,q of the evolution 
through a factor «s(f l o)' f " • 



D. Expansion in ( 

For a two-meson state, the coefficients A n and A' n are functions of the factorization scale u 2 and of the remaining 
kinematical variables C and W 2 . From the definition of GDAs in term of fields given in Eq. (E4j) one obtains moments 



dzz n $ g (z) 



1 



dz 



jn—X 



lp+\n+X 
1 



(P+) 



n+X 



(-^+)"(^(p)^(p')l^)7 + 'Z(0)|0) 
{-id + ) n - 1 (n{ P )n{ P l )\F + ^{x)F+{Q) 



x=0 



x=0 



(53) 



These local matrix elements are the plus-components of tensors that can be decomposed on a basis built up with the 
metric g^ u and the vectors (p + p') M and (p — p')^ . Since (p + p') + = P + and (p — p') + = (2( — 1)P + the moments 
(|S3|) arc then polynomials in 2£ — 1 with degree at most n + 1. The A n and A' n are Gegenbauer moments of ^2 q $ 9 
and <!> ff , respectively, and therefore have the same polynomiality properties in £. Following [ fl9| we expand them on 
the Legendre polynomials, writing 



n+X 

A n ((,W 2 )=6n f J2 B n i(W 2 ) Pi(2C - 1) 

even I 



(54) 



and the analogous expression for A' n with coefficients B' nl . The C invariance properties ( |l7j ) and (26) restrict I to even 
integers in the C even sector. The expansion coefficients B n i arc linear combinations of the local operator matrix 
elements in Eq. (|5^) and are therefore analytic functions in W 2 . As we mentioned in Sect. II C their continuation to 
zero or spacelike W 2 leads to the moments of parton distributions in the pion. 
From Eq. (pi]) the factorization scale dependence of the B n i may be written as 



B nl {W 2 ,f) = B ( +\w 2 ) 



(55) 



with an analogous equation for B' nl involving the factors git . 

In the limit [i — * 00 only the terms with the smallest exponent k[ ^ = in the coefficients (|5l| ) survive. The 
asymptotic form of the distribution amplitudes thus has only n = 1 in the Gegenbauer expansion (|5Cf ) and reads 



E 

9 =1 



$+(z, C, W 2 ) = 18n/ar(l - z){2z - 1) \b{ ] {W 2 ) + b{ 2 \w 2 ) P 2 (2C - 1) 



$ g (z, C, W) = 48z^(l - zf B[ \W 2 ) + B\ 2 \W 2 ) P 2 (2( - 1 



»(-)/ 



(56) 



where P 2 (2C — 1) = 1 — 6£(1 — (). Note that b[ ^ and B\ 2 ' do not depend on a starting scale /j,q because K\ 1 = 0. 

For reasons that will become clear we will also keep the terms with the first nonzero exponent in our model for 

(-) 



(-) 



-(-) 



the GDAs to be developed in Sect. IV. For the quark distribution amplitudes this simply amounts to replacing B 10 
and b[ 2 ' in the first line of Eq. ( jo^ ) with the /i-dependent coefficients Bio & xxd Byi. 

Let us finally remark that, as discussed in |lq| , there is another generalized gluon distribution amplitude, with an 
operator different from the one in Eq. (p4|). It corresponds to pion pairs with angular momentum J z = ±2 and gives 
the leading-twist part of the amplitudes 7*7 — ► tttt where the photon helicities are opposite. The evolution of this 
helicity-two distribution amplitude does not mix with any quark distribution. Its smallest anomalous dimension is 
positive, so that this distribution amplitude tends logarithmically to zero as \x — > 00. The study of this distribution 
would be very interesting. Nothing is, however, known about its size at present, and in our phcnomcnological analysis 
we will neglect its contribution. 
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E. Partial wave expansion 



The decomposition of generalized distribution amplitudes on Legendre polynomials performed in the previous 
section translates into a partial waves decomposition [Hj] if one transforms from polynomials P/(2£ — 1) to p(cos#) 
using that 2( — 1 = [3 cos 8. The rearranged series reads 

nf oo n+l 

5>+ = 6n /Z (l-z)]T Y.Bni(W 2 )C^ 2 \2z-l)P l {cose) (57) 

q=l n=l 1=0 

odd even 

for quarks, where the coefficients B n i(W 2 ) are linear combinations of the form 

B n l = (3 [B n i + Ci, 1+2 B n> 

l-\-2 ~\~ * * • ~r~ Q, n+l B n , n +i] (58) 
with polynomials cj ;/ in 1 . Keeping only n — 1 in the Gegenbauer expansion one is restricted to an S- and a D-wave: 



q=l 



$+ = 18n f z(l-z){2z-l) [B W {W 2 ) + B 12 (W 2 ) P 2 {2Q - 1)] 

= 18n/z(l -z)(2z-l) [Sio(^ 2 ) + Bi 2 (iy 2 )P 2 (cos6l)l (59) 



with 



p 10 (i^ 2 ) = Pio(m^ 2 ) - !^!b 12 (w 2 ), 

Pi 2 (VF 2 ) =/3 2 Pi 2 (W^ 2 ). (60) 

It is a remarkable consequence of the condition I < n+l that the presence of high partial waves implies a departure 
of the two-pion distribution amplitude from its asymptotic form. The ^-distribution of the produced pion pair thus 
contains information about the dependence of the GDAs on z, which as a loop variable is integrated over in the 
amplitude of the process, cf. Eq. (fl2"|). 

One-meson distribution amplitudes are real valued functions due to time reversal invariance. This is not true for 
generalized distribution amplitudes: the two-pion "out" state in the definition (|l^) of $ W7r is transformed into an 
"in" state under time reversal, and these states are different because hadrons interact with each other. Below the 
inelastic threshold, however, two-pion "in" and "out" states with definite angular momentum are related in a simple 
way via the phase shifts of elastic 7T7r-scattering. With the aid of Watson's theorem one then obtains the relation 
B*t = B nl exp(—2i8i) |L9|]. This fixes the phase of the expansion coefficient B n i up to its overall sign: 

B n i = Vm\B n i\ exp(iSi), r\ n \ = ±1, (61) 

where 5i is the tttt phase shift for the l-th partial wave in the 1 = channel. 

F. Momentum sum rule 



Of particular interest are the moments (W,h8|JlS 



1 dz (2z - l)$+(z, C, W 2 ) = (7r+(p)7r-(p0| T++(0) |0), (62) 

II [f ) 

dz<S> g (z,(,W 2 ) = (7r+(p)7r-(p')|r 9 ++ (0)|0), (63) 

where T^ v {x) and T^ v (x) respectively denote the Belinfante improved energy-momentum tensors for quarks of flavor 
q and for gluons. After summing (|6^ ) over all flavors these moments project out the coefficients Biq(W 2 ), Bi 2 (W 2 ) 
and P 10 (W^ 2 ), B; 2 (^ 2 ). 

To proceed one decomposes (7r + (p)7r (p')\ T^ V (Q) |0) on form factors. Their analytical continuation to zero or 
negative W 2 leads to the form factors of the matrix elements (ir + (p)\ Tjf v {Q) \tt + (p')) between one-pion states, with 
W 2 = corresponding to p = p' . At that point we get from Eq. (&32) 
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B 12 (0) = -^-i^, (64) 
9n f 

where is the fraction of light-cone momentum carried by quarks and antiquarks in the pion. No constraint on 
£?io(0) is obtained this way, since the corresponding form factor in the decomposition of (ir + (p) \ T^ v (0) \ir + (p')) is 
multiplied by a tensor that vanishes for p — p'. In an analogous fashion one obtains an expression for -E> 12 (0) from 
the sum rule (|63|). 

We emphasize that both sides of Eq. ( |64| ) depend on the renormalization scale fi. Only the total energy-momentum 
tensor, i.e., the sum T A " y = ^2 q T q ul ' + over quarks and gluons is conserved, so that its matrix elements are 
renormalization scale independent. The appropriate sum of the moments ( |62[ ) and ( p3| ) leads to a linear combination 
of B\2 and B' 12 where the scale dependent term with B[^ indeed drops out and only B\^ ' is left. The normalization 
of (ir + (p) \ T Ml/ (0) \tt + (p)) thus fixes the expansion coefficient 

Bvi\0) = . ( 65 ) 

12 y ' 9n f +48 v ' 

which through the relation (Q) gives the asymptotic value 

< 66 > 

in agreement with the well-known result from the evolution of singlet parton distributions p7fl . 



IV. A SIMPLE MODEL OF THE GDA 



So far no experimental information exists on the two-pion GDA. In the numerical studies to follow we will therefore 
use a simple ansatz for $+(z,£, W 2 ), which is based on the general properties we have discussed in the previous 
section. 

We only consider the contributions from u- and d-quarks, i.e., we take rif = 2. As already mentioned we will use the 
isospin relations ( |2l| ) and (22), and take the asymptotic form of the z dependence given in Eq. (59). It thus remains 
to make an ansatz for the coefficients Bio(W 2 ) and Bi2(W 2 ), or equivalently for Bw(W 2 ) and Bi2(W 2 ) introduced 
in Eq. @. 

For their phases, given by Eq. (^Tj), we use simple parameterizations of the isosinglet S- and D-wave phase shifts 
Sq and S2 obtained in |P8| . They are shown in Fig. ||, where for later use the phase shift S\ of the P-wave is also 
displayed. The result ( pl|) only holds below the inelastic threshold in tttt scattering, therefore we restrict all our 
studies to invariant masses W below 1 GeV. Around that mass, corresponding to the K K threshold, the phase shift 
So of the S'-wave drastically increases. While the analysis of p8fl stops at W — 0.97 GeV and does not exhibit this 
abrupt change, the investigations in Ref. find values of order 200° at W = 1 GeV. Our parameterization of Sq in 
that region is meant to be indicative rather than a precise description of this quantity. Through interference effects, 
the rapid variation of a phase shift leads to a characteristic behavior in the VK-spectrum of appropriate observables 



in our process, as we shall see in Sect. VII 



The analyticity properties of the B n \ and the phase information from Watson's theorem (|61j) may be used to obtain 
the VF 2 -dependence of B n i via dispersion relations, which has been exploited in M,tt9[. Note, however, that while the 
complex phases are simple for the B n i, it is the B n \ that have simple analytic properties in the T4 /2 -plane, given their 
definition through operator matrix elements. The transformation from B n i to B n i introduces extra poles at W 2 = 0, 
cf., e.g., the factors /3 2 — (W 2 — Am^/W 2 in Eq. (|60|). Furthermore, the evaluation of the integrals that solve the 
dispersion relations requires knowledge of the phases at energies above the value of W where B n i is evaluated. This 
further restricts the range of W where B n [ can be obtained using the tttt phase shifts as input. 

To keep our model simple we will make a less sophisticated ansatz. We keep the energy dependent phases and 
S2 from Watson's theorem (|6l|). To determine |-Bio|, I-B12I, and the overall signs 7710, in Eq. (|6l]), we retain only 
the kinematical factors (3 2 in the relation ( |60|) and replace Bio(W 2 ) and -Bi2(W^ 2 ) with their values at W = 0. Close 
to W = 1 GeV one will not expect this to be a good approximation for the S'-wave, given the presence of the /o(980). 
Below this there is however no prominent tttt resonance in the 1 = channel, and the phase shifts show a smooth 
behavior. It seems therefore reasonable to assume that the isosinglet form factors Bio and -B12 do not have a strong 
energy dependence in that region, certainly not as strong as the electromagnetic pion form factor with its large 
variations in modulus and phase due to the p(770). We do however not claim our simple model to be better than, 
say, a factor of 2. 
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FIG. 6. The phase shifts So for the S'-wave, S\ for the P-wave and 82 for the D-wave of elastic tt-k scattering. 80 and 82 refer 
to the 7 = channel. The points are taken from S| , and the curves for 8q and 82 are simple parameterizations. The curve for 
81 corresponds to the parameterization N = 1 of the pion form factor F^(W 2 ) in ]30|]. 
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For the input value of B\2 (0) we use the constraint ( |64| ) with evaluated from the parton distributions in the pion. 
Taking the LO parameterization of GRS |n| we find ranging from 0.5 to 0.6 at a factorization scale between 
1 GeV 2 and 20 GeV 2 . In our numerical studies we use R^ = 0.5. Note that this is very far from the asymptotic 
value (|66|), which for rt/ = 2,3,4 is R^ = 0.27, 0.36 and 0.43, respectively. While using the asymptotic form of the 
^-dependence of the GDA for simplicity (and lack of experimental information) we thus retain a clear non-asymptotic 
effect in the coefficient i?i2(0). We also remark that in the GRS LO parameterization the contribution of strange 
quarks and antiquarks to is at the level of 5% to 10% in a wide range of the factorization scale. This corroborates 
our restriction to u- and d-quarks in the GDA, although with the caveat that the sea quark distribution in the pion 
is not constrained from experimental data |3l) . 

For the coefficient £?io(0) we make use of the relation 



B 10 (0) = -B 12 (0), 



(67) 



which has been obtained in JjiJ using chiral symmetry in the form of a soft-pion theorem. Notice that our ansatz 
then has the property that for (3 — > 1 the S- and D-wave components of the GDA have equal size and opposite sign, 
as is easily seen from Eq. (p0|). 

Putting everything together, we will take the following model GDAs in our numerical studies: 



10z(l- z)(2z- 1)^ 



(68) 



with R v = 0.5. 

With this we can easily calculate the scattering amplitude for 7*7 — * nn to leading order in as. We shall neglect 
here the radiative corrections to the hard scattering, which have been worked out to one loop in |jqj . Taking the 
asymptotic form (^6|) of the quark and gluon GDAs, including the asymptotic value ( |66|) of the ratio R n , they 
were found to reduce the leading-order amplitude for equal photon helicities by 30% if a s — -3, with most of the 



correction being due to the contribution from <f> g . Finally, we recall from the end of Sect. HID that we will neglect 
the contribution of the helicity-two gluon GDA to the photon double helicity-flip amplitude, which is also a one-loop 
effect. 



V. COMPARISON WITH 7*7 -> tt 



Given the close analogy between the production of one and of two pions it is natural to compare the production 
rates of these two processes. Since our estimations for tttt production are at lowest order in as we will compare 
with the corresponding expression for the one-pion case for consistency, although experimental data and more refined 
theory analyses are available there. From the leading-order expression (|l4j) we obtain the cross section for the process 
ej — ► en as 

^ = f^^ (69) 



where we have used the asymptotic distribution amplitude 4>Z = —<^ d — i^/^f^ z(l — z) with f v « 131 MeV. For 



a 



lowest-order approximation, the cross section ( |69| ) is in fair agreement with the data 

To compare with two-pion production, we integrate the cross section for — > e tt^tt 1 from threshold up to W„ 
With our model GDA <M) we find 



da e ^ e ^o _ 25a 3 1 f m< ™ r W 2 L Am l ( \f> 12, l.fi , 
dQ- ~ 72^ QH\~e) J 4ml dW \I 1 'W^{^ +5 |Bl2 ' 



125 a 3 1 2 2 / 4m 2 ( Wl ax 3 mj \ 

- 243 S 2 7 Q 2 (l- e ) Wl ax \±rn% 4 ^J' 1 j 

A consequence of the identical scaling behavior of the two processes is that the ratio of the cross sections (|7(]) and 
( |69| ) is independent of Q 2 in the Born approximation. 

Fig. ^ shows the ratio of the cross sections (|7(]) and ( |69| ) as a function of the upper integration limit W max . We 
see that, even when integrating up to W — 1 GeV, the single- pion production comes out as clearly dominant. We 
remark that the measured production rates M for a single 77 or 77' are comparable to that of a 7r°. With our isospin 
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ej- 
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0.05 - 




W max [GeV] 

FIG. 7. The ratio of the cross sections (^) and ( jpg] ) for the production of 7r°7r° and of ty° in the limit of large Q 2 
cross section for ej —> e7r°7r° is integrated over W from threshold up to W max . 



The 



relation ( pl| ) the cross section for 7*7 — -> 7r + 7r~ is twice that of 7*7 — > 7r°7r°, the relative factor 1/2 for 7r°7r° being 
due to the phase space of identical particles. Due to phase space one does not expect the production of more than 
two pions to be important for W below 1 GeV, except for the decays 77 — ^ 37r and 77' — > hn. The picture thus emerges 
that with our estimation for 7*7 — > 7T7t the production of hadrons in 7*7 collisions up to 1 GeV is dominated by the 
pseudoscalar channel, in other words by the parity-odd sector as opposed to the parity-even one. This is reminiscent 
of the special role played by the axial current in low-energy QCD. 

At this point we wish to comment on the end-point regions of the integrals over z in the factorized expressions 
( fl2| ) and (^) for two-pion and one-pion production. For z — > and z — > 1 the hard-scattering kernels are divergent, 
corresponding to the quark exchanged between the 7 and 7* going on-shell. These poles are canceled by the end-point 
zeroes of the two-pion and one-pion distribution amplitudes, so that the end-point regions give a finite contribution 
to the scattering amplitude in both cases. Quantitatively, the quark virtualities in the hard-scattering diagrams are 
zQ 2 and (1 — z)Q 2 , and it is clear that for a given finite Q 2 there is a region in z where our leading-order expressions 
should receive important corrections. At small virtualities the strong coupling becomes large, increasing the size of as 
corrections, and when zQ 2 or (1 — z)Q 2 becomes comparable to the square of typical transverse quark momenta in a 
pion, then power corrections due to the effect of the transverse momentum of the produced qq- pair will be important. 
We recall in this context that various theoretical attempts to evaluate such corrections lead to fair agreement with 
the data for the 7~7t transition form factor [pp2j] down to rather low values of Q 2 . 

For pion pair production both the hard-scattering kernel and the distribution amplitude are zero at z — 1/2 due 
to the constraints from charge conjugation invariance, so that compared to the single-pion case the integral in z is 
more sensitive to the end-point regions. We thus expect that for intermediate values of Q 2 corrections to the lowest- 
order results will be more important in 7*7 — > 7T7t than they are in 7*7 — > ir°. The experimental comparison of the 
Q 2 -dependence of these two processes will therefore be interesting and may help us to better understand the origin 
of these corrections, which are a subject of considerable importance in the physics of exclusive processes. 

Taking the asymptotic ^-dependence of the distribution amplitudes as an example, we can explicitly see how 
important the end-point contributions are in the leading-order expressions (12) and (|l4|). For single pion production 
the integrand in Eq. ([l4]) is a constant then, so that 50% of the z-integral comes from the regions where z or 1 — z 
is smaller than 0.25. For two-pion production the integrand is proportional to (2z — l) 2 , and 50% of the integrand 
comes from the regions with z or 1 — z smaller than (1 — 2~ 1//3 )/2 w 0.1. Given these numbers, one can expect that 
corrections to our leading-order calculation w ill no t be negligible for Q 2 around 4 GeV 2 , which is the lowest value 
considered in our numerical estimates in Sect. VIII . 



VI. RELATIONS WITH THE PHOTON STRUCTURE FUNCTION 



The exclusive process we consider here contributes of course to the inclusive reaction 7*7 — > X. As we mentioned in 
the previous section, the inclusive process is built up from a limited number of exclusive channels in the mass region 
of W below 1 GeV. Let us examine the connection between our discussion of one- and two-pion production with the 
familiar description of inclusive 7*7 scattering in the kinematical limit we are taking here. 
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The unpolarized cross section for inclusive deep inelastic scattering on a photon, ej — > eX, can be parameterized 
by two photon structure functions Ft and Fl as 



(Z<7, 



C7 — >eX 



2na 2 



1 



dQ 2 dV^ 2 s 2 xQ 2 (l-e 



(2xF T (x,Q 2 ) + eF L {x,Q 2 )) , 



(71) 



where Ft and -Fl respectively give the contribution from transverse and longitudinal polarization of the exchanged 
7*. The transverse structure function Ft is often traded for F2 = 2xFt + Ft- 

At the level of partons inclusive hadron production is described by 7*7 — > qq to leading order in as, which gives 
the well-known expressions [0 



pqq= 3a ^ 



2tt 



In 



l + /3g 
1-/3, 



Q 2 



ijg _ 12a \ ^4 J2 



^^e\x 2 (l-x) 
1 



x 2 + (1 - .t) 2 +4jc(1 - a;)— i 
2m 2 . l + /3 9 



4 



(1 - 2x) 2 + 4x(l 



0, 



W 2 



■hi 



(72) 



where /3 q = (1 — 4m 2 /W 2 ) 1 / 2 . Note that m g is to be understood here as a cutoff parameter, which regulates the 
collincar divergence in the box diagram with massless quarks. 

The limit of large Q 2 at fixed small W 2 we are taking here implies x — > 1 and is different from the Bjorken limit, 
where W 2 is scaled up with Q 2 so that x remains constant. Neglecting terms of order 1 — x ~ W 2 /Q 2 and m 2 q /Q 2 
the expressions (|72|) become 



F qq _ 3a ^ 4 

~ 27T 2- ^ 



In 



1 + & 
1-/3, 



Ff = O 



W 2 

Q 2 



(73) 



We observe that in our limit the leading-order expression for Ft becomes independent of Q 2 , i.e., it has the same 



scaling behavior as the exclusive channels 7*7 



ln[Q 2 /r 



ln[(l 



7r and 7*7 — > 7T7T. This is to be contrasted with the Bjorken limit, 
x)/x] gives rise to the well-known logarithmic scaling violation of 



where ln[(l + f3 q )/(l ~ f3 q )] 
Ft at zeroth order in as- 

Just as in the case of 7*7 — * 7T7t, the contribution Fl from longitudinal photons is power suppressed in our limit. 
Let us add that in the Bjorken limit the hadronic part of Ft, often parameterized using vector meson dominance, is 
only suppressed by a factor In Q 2 with respect to the pointlike part (f72j) , but does not survive our limiting procedure 
here: since hadronic structure functions typically decrease like a power of 1 — x for x — > 1, it becomes a correction in 
W 2 /Q 2 . 

The contribution of our process to the structure functions is, with our ansatz 



for $+, 



T7i7r"'"7r +7r°7r° 
T 



25 a 
96^ 



[3 \B 



10 



625 a 
3456 7T 



Ka 1 



2 9 



(74) 



As a function of W this quickly rises from the threshold at 2m w , levels off for W around 400 to 500 MeV, and then 
remains flat with a value F^, w +7T 17 /a » 0.0077. Let us compare this with the result ( f73| ) of the qq calculation for 
it- and d-quarks (including strange quarks would only lead to a minute change due to the charge factor e q ). With 
the quark masses m u = wid — 290 MeV from the parameterization of the photon structure function by Gordon and 
Storrow Q| we get a value of Fj? /a w 0.15 at W = 1 GeV, much larger than the one we obtain for pion pairs. 

It is worth remembering that 7*7 — > also is the hard-scattering subprocess in our factorized expression for 
7*7 — + 7T7T. As we discussed at the end of the previous section, the collinear divergence of this process shows up as 
singularities at the end-points of the z-integration in Eq. (|l^) and is canceled by the end-point zeroes of the GDA, 
i.e., by the hadronization process. In the calculation of open qq production no such cancellation takes place and 
the divergence of the diagram has to be regulated. This reflects the fact that even in the limit Q 2 — > 00 inclusive 
hadron production from 7*7 cannot be calculated in perturbation theory alone (unlike for instance inclusive hadron 
production from a single timelike photon) and that the separation of Ft into a perturbative pointlike and a non- 
perturbative hadronic part is not unambiguous. While more sophisticated procedures have been elaborated in the 
literature, we consider it sufficient for our purpose to use the quark mass regulator in Eq. (fZ3|) . One might also take 
massless quarks and a lower cutoff k± on the transverse quark momentum, obtaining the same result (|73) with m q 
replaced by kj_ in the expression of f3 q . While keeping us away from the region where perturbation theory breaks 
down, such phenomenological regulators become of course inadequate as one approaches the "threshold" where f3 q = 0. 
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For our numbers this is at W = 580 MeV. One should bear this in mind when using the expression ( |73| ) for invariant 
masses W around 1 GeV. 

On the other hand we saw in Sect. ^ that with our estimate of two-pion production the hadronic mass spectrum 
below 1 GeV is dominated by the single-meson states tt°, rj, i]'. It is clear that in such a region the parton-level result 
can only hold in the sense of parton-hadron duality, averaged over a sufficiently large interval in W. We therefore 
integrate the cross section for — ► eX over W from threshold up to 1 GeV. The parton-level result, obtained from 
Eq. (73) with m u = m e i = 290 MeV, amounts to 2.42 times the cross section ( |69| ) for one-pion production. This factor 
should be compared with the factor 1 + 0.26 + 0.97 + 2.64 for the individual contributions of the exclusive channels 
7r + tttt + ri + if. Here we used Eq. ( |74| ) for two-pion production, whereas for r\ and rf we replaced = 131 MeV 
in Eq. ( |69| ) with the respective values 129 MeV and 213 MeV taken from the analysis of |3^| , Given the caveats of 
parton-hadron duality (below 1 GeV there are very few resonances in the two-photon channel, and W — 1 GeV is 
just above the rf threshold) and those of the parton-level calculation itself (discussed above), we find the agreement 
remarkably fair. 



VII. PHENOMENOLOGY 



We will now discuss the phenomenology of our process in and in e + e~ collisions. The production of neutral 
and charged pion pairs is rather different in this respect, since 7r°7r° is only produced by the 7*7 subprocess we have 
discussed so far, whereas for tt + tt~ production this process interferes with bremsstrahlung, i.e., the production of the 
pion pair from a timelike photon radiated off the beam lepton. We start with the simpler case of neutral pions, and 
then discuss charged pairs. In the following we will restrict ourselves to unpolarized photon and lepton beams. A 
brief discussion of beam polarization will be given in Appendix |b| 



A. Helicity amplitudes 

The building blocks of our investigation are the helicity amplitudes for 7*7 — > tttt, which describe the dynamics of 
this process in a model independent way. They are obtained from the hadronic tensor T^ u by multiplying the reduced 
amplitudes 

Ai 3 (Q 2 ,W 2 ,9) — e^T^e'^, i = + J Q,-, j = +,~ (75) 

with the squared elementary charge e 2 . In the 7*7 cm. our photon polarization vectors read 

£0 = ^(^1,0,0,9°), e± = -^(0 > Tl,-*,0) (76) 



for the virtual and 



e'± = i=(0,-l,+i,0) (77) 



for the real photon, where we have used the coordinate system described in Sect. II A. By parity invariance, there are 

only three independent helicity amplitudes, which we choose to be ^4++, A |_ and Aq + . 

Each of these three amplitudes plays a distinctive dynamical role in the kinematical region Q 2 ^> W 2 ,A 2 . It is 
A ++ that receives the leading twist contribution we have discussed in detail, and which in the scaling regime gives 
access to the generalized quark distribution amplitudes ', 

^ + + = Ef/^^*r(^ 2 ) (78) 

to zeroth order in as- The amplitude A |_ has a leading-twist part at order as, due to the helicity-two gluon GDA. 

We briefly discussed this at the end of Sect. HIE; for more detail we refer to jlfj]. Finally, the contribution Aq+ 



from a longitudinal 7* is nonleading twist. The predicted power behavior in Q 2 at fixed W 2 and C is therefore that 

A++ becomes independent of Q 2 , whereas Aq+ decreases at least like 1/Q- The amplitude A |_ should become Q 2 - 

independent. If the helicity-two gluon GDA is however not sufficiently large, A |_ may be dominated by higher-twist 

contributions at accessible values of Q 2 and should decrease like a power of 1/Q in the corresponding Q 2 -range. 
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Of course all these predictions are to be understood as up to corrections in logQ 2 . At sufficiently large Q 2 , the 

longitudinal amplitude Aq + is thus predicted to be small compared with A ++ . One can also expect that A |_ will be 

smaller than A ++ , since its leading-twist part is suppressed by as- 

To discuss the different partial waves in which the pion pair can be produced, we expand each of the amplitudes 

A +, A_+ as 



Aa= J2 A ljt (Q 2 ,W 2 ) Pji~ l (cos9), 



i = +,0, 



J = 



(79) 



l=j—i 
even 



where P, m denotes the associated Legendre polynomial corresponding to the value of J z of the tttt system in its cm. 



B. The 7*7 subprocess and 7v ir production 



The differential ej cross section for neutral pion pair production reads 



dQ 2 dW 2 d(cos 6>) dip 



1 



a 3 p 
G ~ 16^Q 2 (l-e) 



( 



.4, 



|A__ 



2e |A 04 



cos tp Re {A* ++ A Q+ - A*_ + A 0+ } 2yje(l + e) 
cos2pRe{A* h+ A_ + }2e), 



(80) 

where the subscript G indicates that the pions are produced in a 7*7 subprocess. For 7r°7r° production the phase 
space in Eq. ( |80| ) is understood as restricted to cos9 S (0, 1), ip 6(0, 2ir) because there are two identical particles 
in the final state. We notice the close similarity of the expression (|8(j) with the cross section of the crossed channel 
process of virtual Compton scattering, and much of what we discuss in the following has its counterpart there G3] . 
To obtain the e + e~ cross section we use the equivalent photon approximation 0, 



da„ 



dQ 2 dW 2 d{cos9) dpdx 2 



a 1 

7T X 2 



1+1 



X2) 



In 



r 2 

' max 



(x 2 ) 



V 2 . 



(^2) 



(1-^2) 



dQ 2 dW 2 d{cos6)dp 



(81) 



where Q'^m an d Q'max are the minimal and maximal virtuality of the photon q' , respectively. We have a lower 
kinematical limit Q' 2 nin = x\ m 2 /(X — X2) determin ed by t he electron mass m e , whereas Q'max depends on experimental 
cuts and will be discussed in more detail in Sect. Vfff A . We remark that for a given ee collider energy the variables 
X2 and y are not independent at fixed Q 2 and W 2 , since 



yx 2 



Q 2 + W 2 



(82) 



and that in Eq. ( pl|) one can easily trade dx2 for dy. 

Since the helicity amplitudes Ay are independent of ip they can be partially disentangled from the (^-dependence 
of the cross section, which is completely explicit in Eq. (|o|). In particular, the relative size and the Q 2 -behavior of 
the ^-independent term and of the terms with cost^ and cos2ip allow detailed tests of the scaling predictions. This 
provides indicators on how close one is to the asymptotic regime at finite values of Q 2 . The (^-independent term in 
the large parentheses of Eq. ( |80| ) receives contributions from leading-twist amplitudes and should thus display scaling 
behavior. The coefficient of cos tp is the interference of leading-twist and non-leading twist amplitudes and should be 
suppressed by at least one power of l/Q. Finally, the cos2(/? term should scale or be power suppressed depending on 
the size of the helicity- two gluon GDA. 

Apart from standard fitting techniques a way to separate terms with different angular dependence is the use of 
weighted cross sections. Weighting each event with a function w(ip, 9) we define 



S ei {w) = / dQ 2 dW 2 dfl 



da e j 
dQ 2 dW 2 dfl 



w(p, 9), 



(83) 



where d£l = d(cos9)dp. Notice that since it is not normalized, S e ^(w) is not just the average value of the function 
w(ip, 9), and includes information about the size of the cross section itself. Interpreting S e ~,(w) as a statistical variable 
one can calculate its standard deviation and finds for its relative statistical error (cf., e.g., Q) 
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where 



is the expected number of events for a given integrated luminosity C. Eq. (pi)) generalizes the well-known result that 
the relative statistical error of the cross section, i.e., of S , e7 (l), is X/VN. We emphasize that the method of weighted 
cross sections is very flexible, and that the choice of weights w(<p, 9) can for instance be adapted to experimental 
conditions such as limited angular acceptance or cuts. One can of course take weighting functions that depend on 
other variables than only 9 and tp. In the following we will also use weighted differential cross sections, where only 
some of the kinematical variables have been integrated out while others are held fixed. In a data analysis, one may 
thus use the weighting technique for some variables and fitting for others. 

The weighting technique is convenient to project out different terms in the cross section. As an immediate example 
we note that the terms constant in tp, with cos tp, and with cos 2tp in the cy cross section are obtained from 

dS e7 (cos mtp) 

dQ 2 dW 2 d{cos9) [ ' 

with m = 0, 1 and 2, respectively. If the moments with m = 1 and 2 are measured to be small compared with the 
moment m = 0, this can be because any two of the amplitudes A ++ , A 0+ , A |_ arc much smaller than the third, or it 



may be due to their relative phases. From the theoretical considerations in Sect. VII A the most natural hypothesis 

in this case is however that Aq + and A |_ are small compared with A ++ . 

While the (^-dependence of the cross section ( p0| ) gives access to the various helicity combinations of the real and 
virtual photon, its dependence on 9 contains information on the angular momentum states in which the pion pair 
is produced. A priori there can be arbitrarily high partial waves, but to analyze the ^-distribution in practice one 
will assume that at a give n W only a finite number of them is important, if only for reasons of phase space. It is 
easy to see from Eqs. (|79| ) and ( p0| ) that for a superposition of partial waves I = 0, 2, ... L the moment of cos mtp 
in ( p6| ) is a linear combination of polynomials P^{cos9) with highest degree 2L. Weighting the cross section with 
cos(mtp) P£l +2 (cos9) and integrating over tp and 9 then gives a zero result. Using these weights thus provides one 
way to estimate from experimental data how many partial waves are relevant. Let us recall the physical relevance of 
this information: in the scaling regime the highest partial wave relevant in A ++ provides a constraint on how far the 



two-pion distribution amplitude is from its asymptotic form, as we discussed in Sect. [II E . 

Let us assume that only partial waves with I < L effectively contribute in the cross section ( p0| ) . The ^-dependence 
of the moments in ( j86| ) is then determined by L + 1 coefficients for m = 0, L coefficients for m = 1 and L coefficients 
for m = 2, corresponding to the number of polynomials P^ 1 (cos 9) with I < L. On the other hand, there are 3L/2 + 1 
complex amplitudes Aiji with I < L in the expansion (179), so that there are 3L + 2 real quantities one would like 
to determine. A global phase is however unobservable in the cross section (pCj), and one may for instance refer all 
phases to the phase of A ++ o. The 3L + 1 coefficients one can extract from the dependence of the cross section on 
tp and 9 thus allow one to reconstruct the \Aiji\ and their relative phases. Since the relation between the angular 
coefficients and the amplitudes is quadratic, there will however be multiple solutions in general. More information 
can be obtained with polarized beams, which we briefly discuss in Appendix |b[ 

The situation is simplest if the ^-dependence of the cross section is compatible with the 7r 7r° system being produced 

only in an S- and a _D-wave, and if in addition the (/^-dependence is flat. Assuming that Aq+ and A |_ are negligible 

compared to A ++ , one can then decompose 



Coo + C Q2 P 2 (cos 9) + C 22 [P 2 (cos 9)] 2 (87) 



dQ 2 dW 2 d(cos t 

and project out the coefficients, using that Cw = dS ei (ww ) / (dQ 2 dW 2 ) with weights 

woo = -— (1 - 42 cos 2 9 + 49 cos 4 6), 
16 

35 

W02 = — -5- (1 — 6 cos 2 9 + 5 cos 4 9) , 
8 

35 

W22 = ^ ( 3 - 30 cos2 + 35 cos4 d )- 
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FIG. 8. (a) The coefficients Cw in the differential cross section (^), evaluated with our model GDA ([i8]). They are plotted 
against W instead of W 2 and therefore have been multiplied with a Jacobian 2W. The values of the remaining kinematical 
variables are s ej — 50 GeV 2 and Q 2 = 5 GeV 2 . (b) The same with the alternative ansatz for the GDA described in the text. 



^From these coefficients one can readily extract the amplitudes |A ++ o|, |A ++ 2|, and the cosine of their relative phase. 

In Fig. U (a) we show the coefficients Coo, C32 and C22 obtained with our model GDA (|6|). The interference term 
between the S- and D-waves contains a factor cos(<$o — ^2) and thus is sensitive to the phase shifts. Characteristic 
features in the TV-dependence of C02 are the point where do — 6-2 — 90°, and the sudden change just below TV = 1 GeV 
due to the behavior of the S'-wave. To explore the dependence of these observables on our input GDA we have made an 
ad hoc modification, changing the sign in the prediction (^7|) from chiral dynamics and taking instead -Bio(O) = i?i2(0) 
with i?i2(0) fixed by the constraint (p4[ ) as before. Notice that this flips the overall sign 7710 of the S'-wave in our 
model. The result is shown in Fig. ||(b) and illustrates the sensitivity, especially of the S-D interference, to the 
detailed dynamics of the 7*7 process. 



C. Production of 7r + 7r and interference with bremsstrahlung 



For the production of tt + it~ pairs in e7 collisions, the 7*7 reaction we want to study competes with bremsstrahlung, 
where the pion pair originates from a virtual photon radiated off the lepton cf. Fig. ^. This process produces 
the pion pair in the C-odd channel and hence does not contribute for 7r°7r°. Its amplitude can be fully computed 
for values of IV where the timelike electromagnetic pion form factor F^TV 2 ) is known. The modulus of has been 
well measured in e + e~ — > n + ir~ . By Watson's theorem its phase is equal to the P-wave phase shift Si, provided that 
W is in the range where 7T7t scattering is elastic. This is rather well satisfied for TV up to 1 GeV. In our numerical 
studies we use for the parameterization N = 1 of |30|, which is in good agreement with the data for l-F^I 2 shown 
in Fig. 
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It also gives a fair description of the phase of F^ in the TV-range where we use it, as we see from the 
comparison with the phase shift 8\ in Fig. ^|. 

The contribution of the 7*7 subprocess to the cross section of ej — > e7r + 7r _ has the same form 
ej — > eTT°ir . We recall that with the isospin relation ( |2l| ) the leading- twist helicity amplitude A 
same for neutral and for charged pion pairs. The bremsstrahlung contribution reads 



( |80|) as for 
in ([r|) is the 



da, 



e~f — >e 7T7T 



dQ 2 dW 2 d(cos6) dip 



a 



(3 2(3 2 



B 167T s^ 7 W z e V 



+ cosy? y / 2x(l - x) (2x - l)\/e(l + e) 2 



sin cos ( 



cos 2tp x(l — x) 2e sin 2 9 



(89) 
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Finally, the interference term of the two subprocesses can be written as[] 



do e 



dQ 2 dW 2 d(cos 6>) dtp 



-2ei 



\/2/3 



16tt s 2 ei ^/w 2 Q 2 e{l-e) 



Cq + C\ cos ip + Ci cos 2(p + C3 cos 3ip 



(90) 



with e; = 1 for positrons and —1 for electrons, and coefficients 



C = Rc^F*A ++ } ^2^(1 - xjy/^T+7) cos0 

+Re{f£4o+} (1 - x)^/e{l + e)sme, 
Ci = Re{F;^l ++ } [1 - (1 - a:)(l + e)] sin 61 



Co 



c a 



-Re< 
+Re{F;, 

— Re |f*A_. 
-Re{F;A__ 



V2x(l - x) 2ecos0 
(1 — x) sin0, 
ri/ e(l + e) sin0 



a/2x(1 - ai)\/e(l + e) cos 0, 



(91) 



Remember that in our kincmatical limit 1 — x ~ W^/Q 2 is small. The structure of the bremsstrahlung contribution 
( |S9| ) then becomes rather simple, since at Q 2 3> W 2 the terms in large parentheses reduce to sin 2 0. With the scaling 
predictions for the 7*7 amplitudes discussed in Sect. |VII A| we also obtain the Q 2 -behavior for each of the coefficients 
C n in the interference term (|90|). 

The relative dependence on Q 2 , W 2 and on e of the three contributions to the cross section is controlled by the 
prefactors 



>(l-e)' 



2/3 2 
W 2 



V2/3 



y/W 2 Q 2 e(l - e) 



(92) 



for 7*7, bremsstrahlung and their interference, respectively, and by the pion form factor F^(W 2 ), which appears 
linearly in the interference and squared in the pure bremsstrahlung term. The factors Q 2 and W 2 in ( |92[ ) can be 
traced back to the propagator of the virtual photon in each subprocess, and the extra factor of /3 in the bremsstrahlung 
amplitude reflects the fact that the pion pair is produced in the P-wave there. 

From the factors (^|[) it follows that the 7*7 contribution decreases faster with Q 2 than bremsstrahlung. On 
the other hand the 7*7 process is enhanced at large e, whereas bremsstrahlung profits from small e. To study the 
amplitudes either in the 7*7 contribution or in the interference term, one will therefore go to larger values of e, 
corresponding to small or intermediate values of y (notice that e = 0.8 corresponds to y = 0.5). The behavior in Q 2 
and y of the different contributions to the ee cross section can be seen in Figs, [ll] an d jl2| , respectively, which have 
again been obtained with our model GDA (|(3§|). Notice that apart from the factors (p2|) just discussed, there is a 
global dependence on y and Q 2 through the factor 1/s 2 in the ej cross section and through the variable x-i in the 



real photon flux, cf. Eqs. (|9j) and (82). 

A very strong effect on the relative weight of the different contributions comes from the pion form factor F 7T (W 2 ). 
As one can anticipate from Fig. |l0| it leads to a considerable enhancement of the bremsstrahlung term in a broad 
W interval around the p mass, thereby also enhancing the interference. The VF 2 -dependence of the different terms, 
obtained with of our ansatz (jf^) for the GDA, are shown in Fig. [l3|. As we discussed in Sect. |[y]this ansatz most likely 
oversimplifies the VK 2 -dependence of the coefficients B\q and B12 in <f>+, but the corresponding error in estimating 
the VF 2 -behavior of A ++ should not change the qualitative picture of Fig. ||[ 

In the limit of large Q 2 the different contributions to the cross section have distinctive dependences on ip. The 
7*7 contribution is predicted to be constant in ip with a cos2t/? modulation due to the product , . The 



2 A C program containing the expressions (|80|), (p9[), (po[), (pl[), as well as the amplitude A_| 
GDA (p8f), can be obtained from the authors. 



calculated with our model 
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FIG. 11. (a) The contributions to the differential ee cross section from 7*7, bremsstrahhing and their interference. Kine- 
matical variables are W = 400 MeV, ip = 0, y = 0.1, E x = 3.1 GeV , E 2 = 9 GeV. For the real photon flux in @ we take 
a max _ mrac j an< } _ ^qq MeV as explained in Sect. VIII A. The sign of the interference term corresponds to an e + 7 
subprocess. (b) The same as (a), but with y - 
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FIG. 12. The same as Fig. [H], but as a function of y at Q 2 = 5 GeV 2 
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FIG. 13. The VK-dependence of the different contributions to the differential e + 7 cross section at Q 2 — 5 GeV 2 , ip — and 
s ej — 50 GeV 2 . The corresponding values of y range from 0.1 to 0.12. 

bremsstrahlung term should be flat, and the interference between them should be dominated by cos tp and cos 3<p, 

going with A++ and A |_, respectively. We show examples of the ^-behavior in Fig. remembering that in our 

model A |_ is zero because we have neglected the contribution of the helicity-two gluon GDA. We notice that the 

cos2< / 9 term in bremsstrahlung, which is kinematically suppressed by 1 — x ~ W 2 /Q 2 , is clearly visible at the larger 
energy W — 800 MeV. The ^-dependence, shown in Fig. [l^, is also quite different for the three components of the 
cross section. For the 7*7 term and the interference it depends in detail on the coefficients of the different partial 
waves contributing to the amplitudes A+j. 



The interference between the 7*7 and bremsstrahlung subprocesses provides an opportunity to study the 7*7 
contribution at amplitude level. On one hand this means that one can completely separate the contributions A ++1 

A |_ and Aq + from different photon polarizations. On the other hand it gives access to the phases of these amplitudes 

relative to the phase of the pion form factor F n , which is equal to the 7T7t phase shift 61 in the range of W we are 
considering. In kinematical regions where the bremsstrahlung amplitude is large, especially for W around the p mass 
peak, the interference can also be used to "amplify" the 7*7 signal. 

For this to be useful it is essential that one can cleanly separate the interference term ( pc|) from the pure 7*7 
and bremsstrahlung contributions in the cross section. This is possible since the 7*7 collision produces the pion pair 
in the C-even channel, whereas in bremsstrahlung 7T7t occurs in the C-odd projection. The interference term can 
therefore be separated by reversing the charge of the lepton in the ej collision, a possibility that is automatically 
provided at e + e _ colliders. Alternatively, any observable that is odd under exchange of the ir + and n~ momenta 
is only sensitive to the interference term, which in turn drops out in any observable even under this exchange. In 
terms of the variables we are using, this exchange corresponds to the substitution (9, ip) — > (n — 9, n + ip). This means 
that we have direct access to the interference through the angular distribution of the pion pair in its rest frame. We 
emphasize that on the experimental level this does not require a perfect angular measurement, but only that the 
detection and reconstruction does not introduce a bias between positive and negative pions. 

From the ^-dependence of the cross section one can extract the four coefficients C n in Eq. (pfj|), which determine 

the three quantities Re{F*A ++ }, Re{F*A n+ } and Re{F*A |_}. In fact, they over-determine them, and one can for 

instance use only C%, C2, C3, and keep the information from Co for a cross check. We remark in passing that this is 
owed to the fact that pions have zero spin, otherwise there would be more helicity amplitudes for the 7*7 reaction 
than independent observables one can extract from the ^-dependence. Using the ip- moments Q8q ) with m = 1, 2, 3 
and inverting the relation between C\, C2, C3 and the helicity amplitudes we obtain 



D. Studying the 7*7 subprocess through the interference term 
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FIG. 14. (a) The (^-dependence of the different contributions to the differential e + 7 cross section at Q 2 — 5 GeV 2 , 
W = 400 MeV and y = 0.1. (b) The same for W = 800 MeV. 
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FIG. 15. (a) The ^-dependence of the different contributions to the differential e + 7 cross section at Q 2 = 5 GeV 2 , 
W = 400 MeV, <p = and y = 0.1. (b) The same for W = 800 MeV. 
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K dS ej (w ) 



with a global factor 



and weights 



A <to e7 W f 7T - fl j = 2Re{ Fl A ) sin 2 cos 9, 

y/e(l + e) dQ 2 dW 2 d(cos9) X J I 1 0+ l 

/V " ' /S( ' ^e^n-e} = 2Rc{f;A_+} sinfl (93) 



xe dQ 2 dW 2 d(cos 6) 



, . „,i , _ /" 3 (^) 2 V2 



A^Q^, W V) = -ei — ^— y- , (94) 



2(1 -x) I e fl „„„o.. , 1 - a: ^2 fl , .2. 



w+ = sin 9 cos ip—\ W 2 cos sin cos 2<p -\ (sin # + 4ecos 6) cos 3ip, 

1 " 1 + e xt 



u>o = — sin cos 9 cos 2y> + y y cos # cos 3<y9, 



2(1 -x) /1 + e 



u>_ = — cos3(/?. (95) 

By taking weights that are odd under the exchange of the tt + and tt~ momenta and summing over configurations with 
9 and it — 9 we have canceled the contributions from the pure 7*7 and bremsstrahlung terms in the cross section. We 
remark that our method can easily be adapted to the case where one does not have full acceptance in tp, since the 
moments of cos<^, cos2yj and cos 3^ are always linear combinations of Ke{F* Aij} . 

The functions Wi have been chosen such that they are finite, because the use of unbounded weighting functions is 
problematic. As a consequence, the terms Ke{F*Aij} on the r.h.s. of Eq. ( p^ ) are still multiplied with functions of 9. 
One can avoid the rather strong suppression of angles 9 near or it in Kc{F*A ++ } sin 3 9 if the measurement of the 

moments (^3|) indicates that A |_ and A 0+ are small compared with A ++ . In this case one may replace the weight 

u>+ with cosy, whose moment is dominated by A++sin9 with corrections of order yl — x Ao + and (1 — x)A |_. 

Alternatively, the moment of 



2(1 — x) I e 

w' + = sin 9 cos <p — \ — \ / 2 cos 9 cos 2ip, (96) 

V x V 1 + e 

projects on A ++ sin 2 9 with corrections only of order (1 — x)A |_. In a similar way the moment of cos#cos2c/2 

approximately projects on Aq + sin 9 cos 9 if A |_ is sufficiently small. 

The ^-dependence of the moments (|93| ) contains information on the partial wave decomposition of the pion pair. 
One way to extract the partial waves is of course to fit the ^-dependence of the weighted differential cross sections (|9^) . 
Alternatively, one can use weighted cross sections integrated over both ip and 9. The weight cos 3<p P 2 (cos 9) / sin# 

readily projects out the Zth partial wave in A |_ as we easily see from Eq. (p3|). Note that, since P 2 (cos9) oc sin 2 9, 

this weighting function is a trigonometric polynomial. Similarly, cos 2ip (cos 9)/ sin 9 can be used to obtain the Zth 
partial wave in Aq + if the contribution from A |_ is small enough. 

For A ++ the situation is more complicated, because the functions w + Pi(cos 9)/ sin 3 9, w' + Pi (cos 9) / sm 2 9 and 
cos tp Pi (cos 9)/ sin 9 are all unbounded. The sa me prob lem occurs for the function wq P^ (cos 6')/(sin 2 9 cos 9). In 



practice one may proceed as we discussed in Sect. VII B and restrict the analysis to a finite number of partial waves, 
which has to be determined from the data. Decomposing the coefficient C n in (|9l| ) on polynomials P{ l +l (cos9) one 
can see that if only partial waves with I < L are relevant in the amplitudes Aij , then weighting the cross section with 
cos mp P£ 1 3 (cos 9) and integrating over ip and 9 must give zero. For a restricted number of partial waves one can then 

find weights to project out the corresponding amplitudes. In the case where A |_ and Aq + are negligible and only 

the partial waves I — and I — 2 are important in A ++1 we have for instance 

K_ dS e ~t(w + i 

l-(l-x)(l + e) dQ 2 dW 2 

with 



= Re 



{KA++i}, 1 = 0,2 (97) 
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FIG. 16. (a) Differential cross sections weighted with cos c 
s C7 = 50 G eV 2 and Q 2 = 5 GeV 2 with the model GDA in m 



Sect. VII E 



w+o and w+2. The curves are calculated for an e +r y collision at 
. (b) The same for the alternative GDA described at the end of 



w +0 = — cos ip (1 + 2 cos 

6TT 



16 ,1 A 2 

W + 2 = COS (f (1 — 4 COS 

37T 



(98) 



In Fig. |16| we show the moments of cos tp, w + o an d w +2 as a function of W for our model GDA ( |68[ ) and also for 
the alternative ansatz described at the end of Sect. VII B. We clearly see the sensitivity of our observables to the 
detailed phase structure of the 7*7 amplitude. 



E. Comparison with lepton pair production 

In this section we compare our process e^f — > e n + ir~ with the production of a muon pair, ej — > e , in the same 
kinematics. This is interesting in itself because production is the QED analogue of the reaction we are studying, 

but also because it constitutes an experimental background to the extent that a muon pair can be misidentified as a 
pair of charged pions. 

The helicities of the muons can couple to or ±1 along the direction of the fi + momentum in the 7*7 cm. 
From angular momentum conservation in the subprocesses 7*7 — ► and 7* — > (the latter occurring in 

bremsstrahlung) it is clear that the dependence on 6 and ip must be different in the cross sections for pion and for 
muon pair production. We therefore restrict ourselves here to the cross sections integrated over these angles. For the 
bremsstrahlung contribution we have 



da, 



cj — >eX 



dQ 2 dW 2 

where 



2x(l-x)(l-e) , Xf . 



02 

B ob e~{ 



fs(W z ), (99) 



+7! - f33\FAW 2 )\ 2 +ir 2/3,(3 -fil) 

jb - ^72 . ib - ^72 y im > 

with the muon velocity /3 M = (1 — imZ/W 2 ) 1 ' 2 in the 7*7 cm. For the 7*7 process we can easily adapt the result 
((73|) for open gg-production to the fJ, + n~ case and find 

wl = i^h) fSiw% (101) 
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FIG. 17. (a) The functions f§ occurring in the bremsstrahlung contribution (j9ij) to tt + h and fi + fi production. They 
are plotted against W inste ad of W 2 and therefore have been multiplied with a Jacobian 2W. (b) The same for the functions 
Jq in the 7*7 contribution (101). Note that the curve for pions, obtained with our model GDA (pq), has been multiplied by a 
factor 100. 



where 



G 



V 18 



(3 1 



' G 



8 In 



1 + 



(102) 



up to corrections of order W 2 /Q 2 . 
in the pion and the muon case. 



Notice that both for bremsstrahlung and for 7*7, the Q 2 -dependence is the same 



The functions fg and f@ are compared in Fig. [17[ We see that for the bremsstrahlung contribution pion production 



is enhanced by the strong resonance effect around the p mass, as manifested in F^iW 2 ). In the 7*7 subprocess, on 
the other hand, we find that with our estimate of the GDA, pion production is suppressed compared to muons by a 
factor 50 to 100. This is mostly due to the numerical constants in the expressions (102). In part it also comes from 

the logarithm log(l — /3 M ) in Jq , which is generated by the collinear regions around 9 = and ir as discussed in 



Sect. VI. Notice that for this reason the fx + fjT cross section will be relatively sensitive to cuts that affect 9. The same 
will apply to the interference between bremsstrahlung and 7*7, which drops of course out after angular integration. 
From the results on fg and Jq we expect that the ratio of muon to pion pair production will be appreciable in the 
interference term. 

Another experimental background, again due to particle misidentification, is e ± 7 — > e ± e + e~. Compared with 
production there are further Feynman diagrams, which can be obtained from the muon case by interchanging 
the lines with momenta k' and either p or p' , now corresponding to identical particles. We shall not analyze these 
diagrams here, but will at least assess the contributions from those diagrams that are also present in muon production. 
Repla cing /3 M with e we obtain velocities extremely clos e to 1. Nothing dramatic happens in the bremsstrahlung 
part (100), but the logarithm in the 7*7 subprocess ( |102[ ) is now much larger than for muons. This large logarithm 
is however generated by transverse momenta p± of order m e in the 7*7 cm., which correspond to extremely small 
angles 9 of order m e /W. For any c ut th at effectively leads to a minimum angle 9 m i„ much larger than that, one has 
to replace (3^ with cos 9 m i„ in Eq. (102), which can significantly reduce the size of the logarithm. 

We finally note that the differential cross sections for e + e~ — > e + e~ e + e~ and e + e~ — > e + e~ have been fully 

calculated to first order in QED and are available in the form of Monte Carlo generators [pq] . 
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FIG. 18. (a) The pion energies p° L (solid) and p'® (dotted) in the laboratory as a function of the angle 8 in the two-pion 
cm. The values of the remaining kinematical variables are Ei = 3.1 GeV, E 2 = 9 GeV, Q 2 = 5 GeV 2 , W = 400 MeV, y = 0.1, 
ifi = 0. (b) The same for the transverse pion momenta p ±L and p'± L in the laboratory. 

A. Laboratory kinematics and experimental cuts 

Before giving our estimates for the cross section of our process at various e + e~ colliders, we give a brief discussion 
of the kinematics in the laboratory frame and the effects of some experimental cuts. Starting with the kinematics of 
the scattered lepton k' , we remark that there is a simple transformation between the variables (Q 2 , y) and (E±, an,), 
where E\ and olxl respectively are the energy and scattering angle of k' in the laboratory frame. Imposing minimum 
values on both quantities we have 



1 



AE\ E x 



< 1 



4E? 



j^jfrnin 

Ei 



(103) 



and 



V 



Q 2 1 + cosan 
4E 2 1 — cos a\L 



> 1 



Q 2 1 



cos a 



1L 



4E 2 1 - cos a™ L m ' 



(104) 



The condition (103) cuts on large values of y and is generally not very serious, because mos t inf ormation on the 7*7 

q). The lower cut (104), on the other hand, 
is not large enough. We will encounter an 



process is obtained from low or intermediate y as we discussed after Eq. (Bl 
severely restricts the int eresting y-range in some experimental setups if 
example of this in Sect. VIII B. 

The transformation of the pion momenta into the laboratory system leads to rather lengthy expressions, which we 
will not give here. Notice that the lepton k' has a large transverse momentum k' ±L = Q\J\ — y in the laboratory, 
which must be compensated by the two pions. Even though the tttt system has a rather low invariant mass, the pions 
thus carry large transverse momentum which helps to detect them. An exception are configurations with the the cm. 
angle close to or n, which in the laboratory correspond to an asymmetric sharing of momentum between the two 
pions. This is illustrated in Fig. |lj. 

It is instructive to consider the point where there the tttt system has zero longitudinal momentum in the 
laboratory. With the approximation W 2 <C Q 2 we find 



PL = yEi 



i-y Q 2 

y m' 



(105) 



so that P| 



when y equals 
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For Q <C Ei this simplifies to y — Q/(2Ei). If y is very different from y the tttt system is strongly boosted along 
the beam axis, and if this boost is too large then one or both pions will go out of the detector acceptance. 

We finally have to discuss the kinematics of the scattered lepton I' in the laboratory. In terms of its scattering 
angle ol 2 l we have, up to electron mass corrections, 

l' ±L = (l-x 2 )E 2 sin a 2L (107) 

for the transverse component of V ', and 

Q' 2 = -q' 2 = (1 - x 2 ) E 2 (2sin^i) 2 (108) 
for the photon virtuality. For small a 2 L we obtain the simple relation 

1 — X 2 

It turns out that an antitagging condition on the lepton I', i.e., a 2L < a r 2 "f x with a r 2 "f x determined by the acceptance 
of a lepton in the detector, is not enough to keep Q' 2 small. With the parameters E 2 and a^ x in Tables and III 



we find that, except in the region of x 2 very close to 1, the maximum values of Q' 2 and l' 2 L are a few GeV . Under 
such circumstances it is clearly inappropriate to approximate q' 2 as zero and the momenta I, I' , q' as collinear, which 
we have done throughout this work. Both the kinematical transformation from the ej frame to the laboratory and 
the calculation of the cross section have to be modified then. One must not only recalculate the two-photon and 
bremsstrahlung processes of Fig. |^ but also include further diagrams contributing to the reaction e + e~ — > e + e~ nit. 
Although this is possible in principle, we wish to retain here the simpler expressions for the cross section with one 
real photon. We therefore require that Q' 2 be small compared with the other kinematical invariants in our problem. 



A way to achieve this, suggested by Eq. (109), is to impose an upper cut on l'± L , i.e., in practical terms on the sum 
|kj_ L + p^ L + p' ±L \ of the reconstructed transverse momenta, possibly supplemented by a lower cut on 1 — x 2 . In our 
num erical studies we determine the maximum virtuality Q'^ax in the photon flux of Eq. (|8|) through Eqs. (107) and 



(108) by requiring both a 2L < o^l™ anc ^ l'±.L — 1±l x = 100 MeV. This leads to considerably smaller virtualities than 
the antitagging condition alone, although for x 2 very close to 1 the resulting Q' 2 lax is still not very much smaller than 
W 2 . In practice one may therefore consider an additional cut on x 2 , but we have refrained from this in our estimates. 
Notice that the Q' 2 -spectrum of the photon flux is logarithmic so that a substantial part of the cross section comes 
from Q' 2 much smaller than Q' 2 nax - 



B. B- factories 



We have now all elements to give cross section estimates for existing e + e~ facilities. We start with the -B-factories, 
BABAR, BELLE and CLEO, running at a cm. energy ^/s ee around 10 GeV. Using our model GDA ( |68|) we calculate 
the integrated cross section a and the individual contributions oq and <jb from the 7*7 and bremsstrahlung subpro- 
cesses. To project out their interference term we take simple examples of weighted e + e~ cross sections, 5' ee (sgn(cos if)) 
and 5 ee (cos f), defined in complete analogy with the weighted e^f cross sections (p3[). We remark that S' ee (sgn(cos ip)) 
is simply the left-right asymmetry of the pions in their cm. We integrate over y from its lower kinematical limit 

> Q 2 + W 2 
^ AE\ E 2 

up to y = 0.5. Choosing a larger value increases the cross section, but the gain is mainly due to bremsstrahlung. Up to 
which values of y one can extract useful information on the 7*7 process depends of course on the detailed kinematics 
and must be studied in each particular case. The same is true for the upper limit of the Q 2 -integration. For its lower 
limit we take 4 GeV 2 as a minimum value where one might expect a lowest-order calculation to be reliable, cf. our 
discussion in Se ct. | y|. T o determine the value of Q' 2 iax in the equivalent photon flux we impose the cuts discussed at 
the end of Sect. VIII A. Our results for e + e~ — * e + e~ it + it~ are given in Table |, where apart from the quantities just 



discussed we also give the coefficients in the relative statistical errors 5{w) ~ const/ VN of the weighted cross sections 
S ee (w). We see that the results for the different kinematical situations are practically identical. This indicates that 
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Ei [GeV] 
E 2 [GeV] 
aSf [mrad] 



BABAR 
e~ tagged 
9 
3.1 
684 



BABAR 
e + tagged 
3.1 
9 

300 



BELLE 

e~ tagged 
8 
3.5 
154 



BELLE 
e + tagged 
3.5 
8 
112 



CLEO 

5.3 
5.3 

227 



a [fb] 
°G [fb] 
ob [fb] 

See(sgn(cOSy)) [fb] 
See (COS If) [fb] 

V^V 5(sgn(cos <p)) 
V~N <5(cos (fi) 



452 
15 
437 
-51 
-40 
8.8 
7.7 



452 
15 

438 
52 
40 

8.8 

7.7 



452 
14 

437 
-51 
-40 

8.8 

7.7 



452 
15 

438 
51 
40 

8.8 

7.7 



453 
15 

438 
51 
41 

8.8 

7.7 



TABLE I. Cross sections for e e 
20 GeV 2 , and y from its lower kinema tical lim it (11C) up to 0.5 



real photon flux as described in Sect. VIII A 
and See (cosy) corresponds to a tagged e + 



e+e" n + jr^, integrated over the range W = 300 MeV to 1000 MeV, Q= 4 GeV 2 to 
The cut parameters offi? and l'l£ x = 100 MeV determine the 

sgn(cos ip)) 



In the column for CLEO, the sign of the weighted cross sections S e 



BABAR BABAR BELLE BELLE CLEO 

e~ tagged e + tagged e~ tagged e + tagged 

a?l n [mrad] 300 684 112 154 227 

(tt -a?£ x ) [mrad] 684 300 154 112 227 

[mrad] 300 684 297 524 314 

(tt -0%*°) [mrad] 684 300 524 297 314 

a [fb] 32!) 433 433 443 446 

a G [fb] 6 12 13 13 14 

a B [fb] 323 422 420 430 433 

See(sgn(co S¥ p)) [fb] -31 48 -50 51 52 

See (cosy) [fb] -24 38 -39 40 41 

VN <5(sgn(cos ip)) 10.5 8.9 8.7 8.7 8.6 

\/iV5(cosy) 9.0 7.8 7.6 7.6 7.5 



TABLE II. As Table | but with cuts imposed on the detection angles as specified, and in addition a minimum transverse 
momentum for the tagged lepton and for both pions of 100 MeV in the laboratory. Ei, E2 and ot™£ x for each column are the 
same as in Table fl. 
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it is the cut l' ±L < 100 MeV which determines the real photon flux in most of the relevant parameter space, and not 
the cut on a 2 L , which is different in each of the five cases. We also find that S ee (cos tp) has a slightly smaller relative 
statistical error than S ee (sgn(cos ip)) and thus greater sensitivity to the interference term. 
To estimate the effects of experimental acceptance for the detected particles we impose 

• a cut a™ n < a 1L < ol™£ x on the scattering angle a\L of the tagged lepton k' , 

• a cut (9™ m < (9 L ,6' L ) < Q™ ax on the polar angles 6 L and 0' L of the pion momenta p and p', measured with 
respect to the direction of the initial beam lepton k, 

• a minimum transverse momentum of 100 MeV for the tagged lepton and for each of the pions. 

All quantities refer of course to the laboratory frame. The results are shown in Table [n| 

Comparing with Table | we see that the effects of these cuts are generally quite moderate. The strongest effect 
is observed for BABAR kinematics in the case where t he e ~ is tagged. This can be traced back to the constraint 
a mm <- a xhe minimum value of y implied by Eq. (104) for Q 2 = 4 GeV 2 is 0.46 in this case, which effectively 



cuts away all phase space where the 7*7 process is relevant. The situation improves rapidly as Q goes up, and for 
Q 2 = 6 GeV 2 our cut implies y > 0.19. For the other experimental configurations the same cut is much less restrictive: 
for BABAR kinematics with a tagged e + our cut on ot\L implies y > 0.18 at Q 2 — 4 GeV 2 , whereas in the cases of 



BELLE and CLEO there is not restriction on y from the inequality (104) at all, not even at Q 2 = 4 GeV 2 . 

We find that in the kinematics of B-factories the interference term is clearly larger than the contribution from 
7*7 alone. With several 10 fb _1 integrated luminosity our estimated cross sections give event rates of order 10,000. 
As we see from the tables, the relative statistical error on the interference term, extracted through the moments 
See (sgn(cos if)) or S ee (cosip) is about 8 to 10 times larger than for integrated cross sections (where it is 1/y/N), so 
that the interference could be measured with statistical errors in the 10% range. 

For the production of neutral pion pairs we easily obtain the cross section without cuts by multiplying gq in 
Table | with a factor 1/2, due to the restricted phase space of identical particles. We refrain from a discussion of 
the experimental reconstruction of the four-photon state coming from two pion decays, but for an order-of-magnitudc 
indication of event rates one may take half of the cross sections oq in Table ||. We then estimate hundreds of events 
with several 10 fb _1 , corresponding again to a statistical error around 10%. Thus studies of both charged and neutral 
pair production seem promising to us. 



C. LEP 



Let us now investigate the situation at high-energy colliders, taking as examples LEP1 at E\ — E 2 — 45 GeV and 
LEP2 at E x = E 2 = 95 GeV. 



In the columns labeled "no cuts" in Table IB we list our predicted cross sections, with cuts only on l'^ L and a 2 L so 
that the real photon flux is defined. For the kinematics we have chosen, the cross sections come out about a factor 2 
to 3 larger than at the i3-factories. Luminosities at LEP are however much smaller, so that unfortunately we estimate 
rather low achievable event rates, and it is not clear to what extent studies of our process in this kinematical regime 
will be feasible. 

To see the effect of cuts on the detected particles we require 

• a^l n <a 1L <Tr~ a™l n with a™ L m = 30 mrad and E[ > 0.7 E x for the tagged lepton, 

• 0™ m < (0 L , 0' L ) < 7r — #™ m with 0™ n = 262 mrad, corresponding to pseudorapidities \rj\ < 2, and a minimum 
transverse momentum of 100 MeV for each of the pions. 



The results are given in the columns "with cuts" of Table III. The most serious r estri ction here is the cut on the 
pion angles 9 L and 6' L . This can be understood from our considerations after Eq. ( 106 ). The value of y where the 
7T7T system has zero longitudinal momentum in the laboratory is Q/(2Ei) and thus of order 0.01 to 0.05 here. Over 
most of the j/-range the pions are therefore so strongly boosted in the lab that they appear under extremely small 



angles and cannot be detected. We observe in fact in Table III that the effect of cuts is stronger at LEP2 with its 
higher beam energy, and that it is more pronounced for bremsstrahlung than for the 7*7 process, the latter being less 
affected by a loss of events at larger y. 

At LEP1 the cut on oi\l puts no restriction on y, but for LEP2 we find that for Q 2 = 4 GeV 2 it implies y > 0.5, 
so that one must go to larger Q 2 . For Q 2 of about 8 GeV 2 there is no restriction on y from the constraint (104) any 
more. 
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Ei = E 2 [GeV] 
<?L [GeV ] 


no cuts 
45 
20 


LEP1 


with cuts 
45 
20 


no cuts 
95 
40 


LEP2 


with cuts 
95 
40 


a [fb] 


1023 




167 


1333 




50 


(T G [fb] 


86 




53 


124 




17 


(7B [fb] 


937 




114 


1209 




33 


S e e(sgn(cOS(/s)) [fb] 


128 




41 


159 




13 


S ee (COS^J) [fb] 


101 




32 


125 




10 


VN <5(sgn(cos <p)) 


8.0 




4.0 


8.4 




3.7 


y/N 5(cos ip) 


7.0 




3.5 


7.4 




3.3 



TABLE III. Cross sections for e + e" -> e+e tt+tt , integrated over the range W = 300 MeV to 1000 MeV, Q = 4 GeV to 
Qmax, an d y from its lower kinematical limit (11C) up to 0.5. The col umns marked "no cuts" correspond to imposing only the 



cuts that determine the real photon flux as explained in Sect. VIII A, with parameters a™£ x = 30 mrad and 1™£ X = 100 MeV. 
The columns "with cuts" refer to the additional cuts described in the text. The sign of the weighted cross sections S ee (sgn(cos ip)) 
and See (cos tp) is for a tagged e + . 



We finally note that at the very large values of Q 2 accessible at high-energy colliders one can afford invariant masses 
W well above 1 GeV, while still fulfilling the basic condition W 2 <C Q 2 of our study. We have not explored this mass 
region, since our model for the pion GDA is not applicable there. It is however clear that there will be a strong 
enhancement of the GDAs at W around the masses of C-even resonances, such as the /2(1270). 



IX. SUMMARY AND OUTLOOK 



In this paper we have analyzed in detail the process 7*7 — * 7T7t in the domain where the virtuality Q of the 7* is much 
larger than the invariant mass W of the two-pion system. It factorizes into a parton-level subprocess, which is under 
perturbative control, and non-perturbative matrix elements called generalized distribution amplitudes. This makes 
the reaction a laboratory to study the non-perturbative dynamics of a two-pion system forming from a well-defined 
partonic state, namely from a quark- antiquark or a two-gluon pair produced at small distance. The perturbative stage 
of the overall process is completely analogous to the one in single- meson production, well studied in the case of a 7r°, 
T) and rf . It results in a scaling behavior of the amplitude as Q 2 increases at fixed W 2 , selects characteristic helicity 
combinations of the two photons, and predicts that the two pions are produced with total isospin zero. The dynamical 
content of the non-perturbative matrix elements, on the other hand, is more complex than for a single particle. Even 
the lowest Fock state of \ir) (g> that is, qq ® qq, contains more partons that the initial qq or gg system from which 
the two pions are formed. In this sense a GDA describes the transition between different parton configurations in 
the non-perturbative regime. The two-pion distribution amplitude contains the full strong interactions between the 
two pions, leading to dynamical phases which, by Watson's theorem, are identical to the phase shifts in elastic irir 
scattering as long as W is below the inelastic threshold. We use this relation as an input for our model GDA, and 
therefore restrict our study to the W-region up to 1 GeV. 

The evolution equation giving the factorization scale dependence of the GDAs is more complex than for a single pion 
due to the mixing of qq or gg amplitudes, and we have given the relevant splitting functions and anomalous dimensions 
for the quantum numbers of relevance here. A simultaneous expansion of $(2, C, W 2 ) in the parton momentum fraction 
z and partial waves of the pion system leads to local matrix elements between the vacuum and a two-pion state. By 
analytic continuation they are related to the moments of the parton distribution functions of the pion. We have used 
the quark momentum fraction in the pion, determined from a global fit of these distributions, as an input for our 
model of §(z, £, W 2 ). The corresponding value of is well below its asymptotic value under perturbative evolution, 
which may be an indication that the lowest non-asymptotic terms in the crossed-channel quantity ^(z, C, W 2 ) are not 
small at factorization scales in the GeV range. We emphasize that the question of how close one is to the asymptotic 
result of evolution is particularly interesting, because in the case of light pseudoscalars the single-meson distribution 
amplitudes may be surprisingly close to their asymptotic form even at low scales J^j32|. 

From a theory point of view it is also interesting to consider & q (z, (, W 2 ), defined by the matrix element in Eq. fllq), 
for values of W much larger than the scale of non-perturbative interactions. While the dynamics in <E> g (z,£, is 
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entirely soft for small W, part of it becomes hard when W increases. In the limit If >1 GeV and to leading order in 
as one can explicitly write Qq(z, £, W 2 ) in terms of a perturbative subprocess and the qq distribution amplitudes for 
each separate pion [B7|. The resulting Q q (z,(, W 2 ) is very far from the asymptotic form in z. It receives substantial 
contributions from high partial waves of the tttt system, has a power- law falloff like 1/W 2 , and its imaginary part is 
small compared to its real part. 

We have constructed a model for the GDA at W below 1 GeV, using simple structure as a guide, and and 
the tttt phase shifts as phenomenological inputs. Comparing the rates for the production of tttt and of a single 
pseudoscalar meson, we found that the hadron spectrum in 7*7 collisions below 1 GeV is strongly dominated by the 
single resonances 7r°, 77, and rf . 

We have further compared our process with open qq production, which at higher invariant masses W is commonly 
used to describe the part of the total hadronic 7*7 cross section due to the pointlike part of the real photon. Interest- 
ingly, we find that in our particular kinematical limit, the corresponding scattering amplitude has the same scaling 
behavior and helicity structure as the one for the exclusive processes 7*7 — > tt and 7*7 — > irir. The main difference is 
that in the tt and tttt cases the collinear divergence of the lowest-order hard scattering diagrams is regulated by the 
hadronization process. This is encapsulated in the distribution amplitudes, which vanish at the end points z = and 
1. In the open qq calculation, on the other hand, the divergence has to be regulated explicitly. We also note that the 
sensitivity to the soft end-point region may be larger for pion-pair production than for a single pion, because for two 
pions the hard scattering and the distribution amplitudes vanish at z = 1/2 for symmetry reasons. Thus one may 
expect the onset of the scaling behavior to occur at different Q 2 in the two cases, an issue that will be interesting to 
study in experiment. 

An investigation of the structure of the cross section shows that in cy and e + e~ collisions information on the 
7*7 process can be obtained either through the square of the 7*7 amplitude, or from its interference with the 
bremsstrahlung process if the pions are charged. This interference can readily be projected out by appropriate C-odd 
observables, and it offers the opportunity to separate the different 7*7 helicity amplitudes. It further provides direct 
access to their dynamical phases, although a full phase reconstruction requires polarized beams (cf. Appendix^). 

The angular distribution of the pion pair in its cm. contains detailed information about the dynamics of the 7*7 
process. The dependence on the azimuth tp separates the different helicity combinations of the real and virtual 
photon, each of which plays a distinct role in the scaling limit. In particular it permits one to study leading-twist and 
non-leading twist amplitudes at the same time, which should provide additional insight into how far one is from the 
asymptotic regime. The ^-dependence, on the other hand, gives access to the partial waves in which the two pions are 
produced. It is sensitive to the phases, which reflect the dynamics of the tttt system and its resonances. Even though 
one will probably not be able to perform a full extraction of the tttt phase shifts in this way, our process provides 
constraints on these quantities that are independent of the analyses of elastic tttt scattering. The presence of higher 
partial waves would in itself be very interesting, since it gives indirect information on the deviation of &(z, £, W 2 ) 
from its asymptotic form in z. 

We have restricted ourselves to the production of pion pairs in this work, but it is clear that many of our results 
are also valid for other exclusive systems. The most obvious generalization is to charged or neutral KK pairs, whose 
comparison with tttt would allow one to study aspects of flavor ST/ (3) breaking in the context of the quark-hadron 
transition. At even higher values of W 2 there is the production of pp, where extra spin degrees of freedom come in, 
as in the well-studied case of the parton distributions of the nucleon. 

Another very similar process is the production of p + \i~ pairs, i.e., the QED analogue of our reaction. Comparing 
the rates of e"f — > ep + p~ with our estimate for ej — > eTT + TT~ we find that the bremsstrahlung mechanism prefers 
pions if W is in the vicinity of the p mass, reflecting the strong resonance effect in the tttt system. For the production 
from 7*7, on the other hand, the cross section is considerably larger in the case of muon pairs. We remark that this 
could not be anticipated from a dimensional analysis. The amplitudes for 7*7 — > pp and for 7*7 — > tttt have the same 
Q 2 -dependence in our kinematical limit, and the two-pion distribution amplitude, which describes that pions are not 
pointlike but have internal structure, is a dimensionless quantity. 

Using our model GDA to calculate the cross section for e + e~ — * e + e _ tttt, we find encouraging rates for the 
kinematics and luminosity of i3-factories. Thus there should be enough statistics for detailed studies at these facilities. 
Our estimates of the effect of cuts also indicate that in the kinematical region interesting in our context, the pions and 
the tagged lepton are well within the experimental acceptance. For high-energy colliders such as LEP, our predictions 
are less optimistic, at least in the range of W below 1 GeV which we have studied here, due both to the lesser 
luminosity and the strong longitudinal boost of the pion system. 

In conclusion, we find that the process 7*7 — > tttt can offer valuable insight into the interactions between quarks, 
gluons and hadrons, and that it should well be measurable at existing e + e~ facilities. 
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APPENDIX A: PION ISOSPIN STATES 

We specify in this appendix our sign convention for the definition of pion states. The relative sign for 7r + and ir~ 
is relevant because it determines the relative sign of the GDAs for charged and neutral pion pairs. 
In terms of eigenstates \ir l ) of the isospin operators P (i — 1, 2, 3) we define 

K+) = -L (1*1) + , = ± (i^i) _ i)^)) , k o } = ^ (A1) 

Notice that the sign for \tt + ) is opposite to the usual convention for eigenstates of SU(2). This has to be remembered 
when writing down two-pion states with definite isospin using the Clebsch-Gordan coefficients. 

The convention ( |Al[ ) is in line with the customs of field theory, see for instance Sect. 12.5 of |38| ]. If, starting from 
the real scalar fields associated with {ir 1 ) and \n 2 ), one constructs the complex scalar field <p which creates |vr ) out of 
the vacuum, then is created by the conjugated field ip*. If one used the opposite sign in defining |7r + ), which is 
more natural in the context of isospin, then there would be an extra minus sign between the fields creating 1 7r ) and 



\tt + ). Through the LSZ reduction formula this sign would show up in crossing relations. With our definition (Al) 
this does not happen, and we have for instance that the spacelike pion form factor 

{n + (p)\JU0)\* + (p')) = (p + p'TFAt) (A2) 

with t = (p — p 1 ) 2 becomes 

(n+(p)n-(pi)\ JUO) |0> - {p-pYFM (A3) 

with s = (p + p') 2 in the timelike region. We remark in passing that if one uses the isospin relation (|2^ ) and neglects 
the contributions from strange and heavy quarks, one has the sum rule 

dz<i>-(zX,W 2 ) = (2(-l)F 7r (W 2 ). (A4) 



The choice ( Al) also leads to a convenient relation for the action of the charge conjugation operator C, namely 

C\n+) = |7r->, C|tt°) = |tt ). (A5) 

The impossibility to find a sign convention that is natural for both charge conjugation and the isospin algebra is 
discussed at le ngth in Chapt. 5, §7 of pjJ (where the other sign in defining |7r + ) was chosen). We also remark that 



the definition (Al) implies 

(7r + |« a (x)d^(0)|0) = <7r-|4(a?)«/j(0)|0), (A6) 

and therefore a relative plus sign between the distribution amplitudes for 7r + and 7r~ . 

Our definition is the same as the one chosen by Polyakov et al., cf. |lq |, and it was also adopted in j37j. We finally 
mention that the definition leading to Eq. (15) of Q] has the opposite sign for \tt + ). 

APPENDIX B: BEAM POLARIZATION 



As we have shown in Sects. VII B and VII D, the unpolarized e'y cross section contains detailed information on the 
7*7 helicity amplitudes A^j. From Eqs. (|80|) and ( |9l| ) it is however clear that this information is not sufficient to 

fully reconstruct the three independent complex amplitudes A ++ , Aq + and For completeness we give in this 

appendix the expressions of the cross section with longitudinally polarized lepton and photon beams, and discuss 
what additional information can be obtained from single and double polarization asymmetries. 



3G 



Starting with the 7*7 contribution, we have 



da, 



e*y — >e tttt 



dQ 2 dW 2 d{cos9)dip 



G 



° q - (80) + T&tL Q 2 {l-e) ( Pl Sin ^ Im ( A *++ A °+ A -+ A o+) 2 V^0^) 



P 7 sin tp Im {A* ++ A 0+ + A*_ + A Q+ } 2y^(TTe) 
P 7 sin2<^ lm{A* ++ A_ + } 2e 



PP 7 {i^+i 2 -!^! 2 }^!^ 



-P/P 7 cos ^ Re {^;+^ 0+ + A*_ + A 0+ } 2^e(l - e) J , (Bl) 

where p and P 7 respectively denote the longitudinal polarization of the lepton and photon beam, ranging from — 1 
to 1. Together with Eq. (|8^) we see that if both lepton and photon are polarized, one has enough independent terms 

to reconstruct the real and imaginary parts of the interferences |_, A* ++ A 0+ and A*_ + A Q+ . Furthermore, the 

squared terms |A ++ | 2 and \A |_| 2 come with a different relative sign in the unpolarized cross section and the double 

polarization asymmetry. 

The bremsstrahlung contribution to the cross section reads 



da. 



cy — >e tttt 



dQ 2 dW 2 d{cos6) dip 



eq. (89) + — 



16tt s 2 W 2 e 

L 1 



|P.(VK 2 )| 2 PP 7 ( (2x-l)y/ 



1 - e 2 sin 2 1 



+ COS tp \/2x(l - x)y/e{l - e) 2 sin cos J . (B2) 



Notice that it only contributes to the unpolarized cross section and the double polarization asymmetry, but not to 
single polarization asymmetries. Finally, the interference term can be written as 



da* 



dQ 2 dW 2 rf(cos 0) dp 



eq. (90) - 2e; ^- 4" 



V2P 



16tt s 2 7 y/W 2 Q 2 e{l-e) 



Pi (c[ sin p + C\ sin 2p 
P 7 ( C*7 sin p + sin 2p> + C3 sin 3<p 



Pi P 7 ( + C*J 7 cos p + C2 7 cos 2<p 



(B3) 



with coefficients 



C[ = -Im|F;A ++ |xv / l - e 2 sin0 
+Im{p;A_+} (1 - x)Vl-e 2 
Ci = lm{F;Ao+}xy/e(l-e) 



sm t 



-Im|p;^_+} ^V^O - e)cos0 



(B4) 



for lepton polarization, 



-Im 



C2 



c 7 



{p;A ++ }[l-(l-x)(l-e)]sin0 
+Im|p*A 0+ | ^2^(1 - x) 2ecos0 
{p;A_ + } (l-s)Binfl 
{p;A o +}x v /e(l + e)sin0 
+Im|p;,4_ + | V^(l - ^-^(l + e)cos0, 

im{p;A_ + } 



— Im 
Im 



xe sm f 



(B5) 
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for photon polarization, and 



C l J = Re{i^A ++ } y/2x{l-x)y/e{l-t) cosO 

-Re{i^ 0+ } (1 - x)y/e(l-e) sin 9, 
C[ 7 = Re{F;A ++ } xy/l -e 2 sinfl 

-Rc{f;A_+} {l-x) v / l-e 2 sine, 
C2 7 = -Re jf* Ao + 1 a;-\/ e(l - e) sin0 

+Rc|f; A_+ 1 y/2x(l - x) v/e(l - e) cos 6 (B6) 

if both lepton and photon are polarized. We see that with polarized photons one can extract Im{F^^4 ++ }, lm{F*Ao + } 

and Iva{F*A [_}, which together with the unpolarized interference term makes it possible to reconstruct all three 

complex 7*7 amplitudes for values of W where the pion form factor F v is known. One cannot achieve the same with 
a polarized lepton beam alone, since there are only two terms in the ^-dependence. In this case one can still use the 
suppression by 1 — x of the second term in C\ in order to approximately extract lm{F*A ++ }. Finally, the double 
polarization asymmetry gives access to the same quantities one can already obtain in the unpolarized case. 
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